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Abstract. We consider the anisotropic Calderon problem of recovering 
a conductivity matrix or a Riemannian metric from electrical boundary 
measurements in three and higher dimensions. In the earlier work [S], 
it was shown that a metric in a fixed conformal class is uniquely deter- 
mined by boundary measurements under two conditions: (1) the metric 
is conformally transversally anisotropic (CTA), and (2) the transversal 
manifold is simple. In this paper we will consider geometries satisfying 
(1) but not (2). The first main result states that the boundary mea- 
surements uniquely determine a mixed Fourier transform / attenuated 
geodesic ray transform (or integral against a more general semiclassical 
limit measure) of an unknown coefficient. In particular, one obtains 
uniqueness results whenever the geodesic ray transform on the transver- 
sal manifold is injective. The second result shows that the boundary 
measurements in an infinite cylinder uniquely determine the transversal 
metric. The first result is proved by using complex geometrical optics 
solutions involving Gaussian beam quasimodes, and the second result 
follows from a connection between the Calderon problem and Gel'fand's 
inverse problem for the wave equation. 
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1. Introduction 

The anisotropic Calderon problem consists in determining the electrical 
conductivity matrix of a medium, up to a change of coordinates, from cur- 
rent and voltage measurements made at the boundary. More generally the 
problem may be posed on a smooth Riemannian manifold with boundary. In 
this case the question is to determine the geometric structure of the manifold 
from the Cauchy data of harmonic functions. The purpose of this paper is 
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to study the anisotropic Calderon problem in transversally anisotropic ge- 
ometries, where the manifold admits a distinguished Euclidean direction, 
and to prove uniqueness results for inverse problems in this setting. 

Let {M,g) be a compact oriented Riemannian manifold with smooth 
boundary dM. Harmonic functions in M are solutions of the Laplace- 
Beltrami equation 

AgU = in M. 

Here, the Laplace-Beltrami operator is given in local coordinates by 

OXj \ OXj J 

where (gjk) is the metric in local coordinates, {g^'') = {gjk)~^, and \g\ = 
det{gjk)- Here and below we are using the Einstein summation convention. 

The boundary data of harmonic functions on M is given by the Cauchy 
data set 

Cg = {{u\9M,d^u\9M) ■,AgU = in M, u G H\M)}. 

The normal derivative di^uloM = {du,i')\QM, where v is the 1-form corre- 
sponding to the unit outer normal of dM, is interpreted in the weak sense 
as an element of H'^/'^idM). It is clear that if i/j : M M is a diffeomor- 
phism satisfying ^p\^M = Id, then C^*g = Cg. On manifolds of dimension 
> 3, the anisotropic Calderon problem [26] amounts to proving that Cg 
uniquely determines g up to isometry. 

Conjecture. Let (M,gi) and {M,g2) be two compact Riemannian mani- 
folds with smooth boundary, and let dim(M) > 3. If Cg-^ = Cg^, then 

92 = 1/^*91 

where -0 : M — )• M is a diffeomorphism with V'laA/ = Id. 

This statement has only been proved for real-analytic metrics |26j with 
topological assumptions relaxed in [21] , [2S] , and for Einstein metrics (which 
are real-analytic in the interior) [T^]. The general case remains a major open 
problem, and we refer to [8] for a discussion and further references. The 
corresponding two-dimensional result, involving an additional obstruction 
arising from the conformal invariance of the Laplace-Beltrami operator, is 
known [25] . See [3], [5] for another interesting approach to this problem. 

The work [8] introduced methods for studying the anisotropic Calderon 
problem in manifolds which are not real-analytic, but where the metric has 
certain form. This was based on the concept of limiting Carleman weights, 
introduced earlier in the Euclidean case in [23]. One of the main results of [8] 
states that on a simply connected open manifold, the existence of a limiting 
Carleman weight is equivalent to the existence of a nontrivial parallel vector 
field for some conformal metric. Locally, this condition is equivalent with 
the manifold being conformal to a product of a Euclidean interval and some 
(n — l)-dimensional manifold. We formalize this notion in two definitions: 

Definition. Let (M, g) be a compact oriented manifold with C°° boundary. 
In this paper we always assume that n = dim(M) > 3. 



CALDERON PROBLEM IN TRANSVERSALLY ANISOTROPIC GEOMETRIES 3 

(a) {M,g) is called transversally anisotropic if {M,g) CC {T,g) where 
T = M X Mo, g = e® go, (M, e) is the Euclidean line, and (Mo,g'o) 
is some compact (n — l)-dimensional manifold with boundary called 
the transversal manifold. 

(b) (M, g) is called conformally transversally anisotropic (CTA) if (M, eg) 
is transversally anisotropic for some smooth positive function c. 

Examples of CTA manifolds include compact subdomains of the model 
spaces W^, sphere 5"" minus a point, or hyperbolic space compact sub- 
domains of locally conformally flat manifolds such as 3D symmetric spaces 
as long as they are contained in a conformally flat coordinate neighborhood, 
and conformally warped products 

{M,g) CC{Rx Mo, g), g = c{e®fgo) 

where / is a positive function depending only on the Euclidean variable in 
M X Mq. See tSj, [Hj, t28j for more details. 

The first main theorem in this paper considers the anisotropic Calderon 
problem in a fixed conformal class. Since any conformal diffeomorphism 
fixing the boundary must be the identity map, there is no obstruction to 
uniqueness arising from isometries in this case (see [29] ) . The article [8] gave 
a uniqueness result for this problem on CTA manifolds if additionally the 
transversal manifold {Mo, go) is simple, meaning that any two points in Mo 
are connected by a unique geodesic depending smoothly on the endpoints 
and that dMo is strictly convex (its second fundamental form is positive 
definite). Moreover, a reconstruction procedure was given in [22]. The 
proof used the fact that the geodesic ray transform is injective on simple 
manifolds. On general transversal manifolds we use the following definition. 

Definition. We say that the (geodesic) ray transform on the transversal 
manifold {Mo, go) is injective if any function / G C{Mo) which integrates to 
zero over all nontangential geodesies in Mq must satisfy / = 0. Here, a unit 
speed geodesic segment 7 : [0,1/] — >• Mq is called nontangential if 7(0), '^{L) 
are nontangential vectors on dMo and ^{t) G M™* for < t < L. 

Theorem 1.1. Let {M,gi) and (M, 172) be two CTA manifolds in the same 
conformal class. Assume in addition that the ray transform in the transver- 
sal manifold is injective. If Cg-^ = Cg^, then gi =92- 

In fact this result is a consequence of a corresponding result for the 
Schrodinger equation. Let q £ L°^{M), and define the Cauchy data set 
for the Schrodinger operator — + 9 by 

Cg,q = {{u\aM,duu\aM) ; (-A,, + q)u = in M, u£ H\M)}. 

Again, the normal derivative di^u\dM is interpreted in the weak sense as an 
element of H-^/'^{dM). 

Theorem 1.2. Let {M,g) be a CTA manifold, and let qi,q2 G C(M). 
Assume in addition that the ray transform in the transversal manifold is 
injective. If Cg^g-^ = Cg^g^, then qi = q2- 
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Starting from the pioneering works [7], |41j . the standard approach for 
proving uniqueness results for the Calderon problem is based on special 
complex geometrical optics solutions to elliptic equations. The paper [8] 
presented a construction of such solutions on CTA manifolds and proved 
Theorems 11.11 and 11.21 under the additional restriction that the transversal 
manifold (Mo,g'o) is simple. Some parts of the argument, for instance the 
Carleman estimates required for the construction of correction terms, were 
valid without this additional restriction. However, in the end the simplicity 
assumption was used to produce solutions that concentrate near geodesies in 
(Mo, go) and also to show that the potentials can be determined by inverting 
the geodesic ray transform (actually with attenuation) in the transversal 
manifold. 

In this paper we remove the simplicity assumption on the transversal 
manifold in the construction of complex geometrical optics solutions, and 
prove Theorems 11.11 and 11.21 on any CTA manifold for which the ray trans- 
form is injective. Injectivity of the ray transform is known to hold in the 
following classes of manifolds {MQ,go): 

(a) Simple manifolds of any dimension (see [37j). 

(b) Manifolds of dimension > 3 that have strictly convex boundary and 
are globally foliated by strictly convex hypersurfaces ([44J). 

(c) A class of non-simple manifolds of any dimension such that there are 
sufficiently many geodesies without conjugate points and the metric 
is close to a real-analytic one (see [39j for the precise description of 
this class). 

(d) Any manifold having a dense subset that is covered by totally geo- 
desic submanifolds in which the ray transform is injective (injectiv- 
ity of the ray transform follows immediately from the injectivity in 
the totally geodesic submanifolds). Examples include subdomains 
of (A'^i X N2,hi(B 112) where (A'^i, hi) has injective ray transform and 
{N2,h2) is any manifold. 

(e) There are counterexamples to injectivity of the ray transform. The 
standard one is the sphere with a small cap removed: any function on 
the sphere that is odd with respect to the antipodal map and vanishes 
near the removed cap integrates to zero over nontangential geodesies. 
See also [2] , [JO] for microfocal analysis of the ray transform in some 
non-simple geometries. 

In fact. Theorems 11.11 and 11.21 involving the ray transform will be ob- 
tained as a special case from a more general complex geometrical optics 
construction on CTA manifolds. If {M,g) is a CTA manifold, so {M,g) CC 
(M X MQ,g) for some compact manifold (Mo,(7o) where g = c(e © go), we 
denote points on M by x = {xi, x') where xi is the Euclidean variable and 
x' G Mq. If g G L°^{M), we will consider solutions of the Schrodinger 
equation (— + q)u = in M of the form 

u = e''''^c~^{vs + rs) 
where s is a slightly complex large frequency, 

s = T + iX, 



CALDERON PROBLEM IN TRANSVERSALLY ANISOTROPIC GEOMETRIES 5 



where the real parameter r wih tend to infinity while A G C is fixed, and 
where Vg = Vs{x') £ C^(Mo) is a quasimode, or approximate eigenf unction, 
with frequency s in the transversal manifold. The correction term will 
satisfy ||'^s||L2(Af) — )■ as r — oo. The concentration properties of the 
quasimodes Vg in the high frequency limit as r — )■ oo will be crucial in 
determining properties of the potential. 

Definition. If A G C, we denote by A^a the set of all bounded measures /i 
on Mo for which there is a sequence {Tj)Jl-^ with tj — )• oo and a sequence 
C C^{Mo) satisfying 

II (-^90 - (Tj + i>^f)vj\\LHMo) = 0{Tj), \\vj\\L2(Mo) = 0(1) 

as j — )• OO, such that in the weak topology of measures on Mq one has 

lim \vj\'^ dVgg = n 

where dVg^^ is the volume form of (Mo,go)- 



Theorem 1.3. Let {M,g) be a CTA manifold, and let qi,q2 € C{M). If 

d^i{x') = 



^9,91 ~ ^9,92! then 



oo 

e-^i^-^{c{qi-q2)){xi,x')dxi 



— oo 



I Mo 

for any A G C and any ^ ^ Mx. Here qi — q2 is extended by zero to M x Mq. 

The measures fi £ Aix are called semiclassical defect measures, or quan- 
tum limits, of the families of quasimodes {vr+ix)- The properties of such 
measures are the central object of interest in the study of high frequency 
limits of eigenfunctions and in quantum ergodicity. In general, the dynamics 
of the geodesic flow of the underlying manifold (Mo, (70) will be visible in the 
semiclassical measures. These topics have a large literature, and we refer 
to [T7], [IS], [l6] for surveys. However, our situation seems to be somewhat 
different from most of these works for the following three reasons: 

(1) We only have access to limit measures in the base manifold Mo 
instead of the more usual phase space measures in T*Mo. 

(2) The measures Mx are associated to a family of quasimodes in a 
manifold (Mo,go) with boundary, but there is no boundary condi- 
tion imposed on the quasimodes. This leads to a certain amount of 
flexibility in our setting. 

(3) It is useful to consider measures for slightly complex frequencies 
r + iA where Re(A) is nonzero. 

Theorem 11.21 will be obtained from Theorem 11.31 by a rather direct con- 
struction of Gaussian beam quasimodes that concentrate on a given non- 
tangential geodesic. This construction goes back at least to [1], [11] and has 
been developed by many authors (often for hyperbolic equations), see for 
instance [18], [32]. In our case, the relevant result is as follows. The fact 
that the frequency is slightly complex leads to the attenuated geodesic ray 
transform with constant attenuation — 2A, but eventually analyticity will 
allow to make a reduction to the case A = 0. 
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Theorem 1.4. Let (Mo,go) be a compact oriented manifold with smooth 
boundary, let 7 : [0, L] — )• Mq be a nontangential geodesic, and let A G M. For 
any X > there is a family of functions {vs) C C°°(Mo), where s = t + iX 
and T > 1, such that 

||(-Ago - s^)t;,||i2(Mo) = 0(t-^), |bs||L2(Mo) = 0(1) 
as r — )• CO, and for any ^ S C(Mo) one has 

lim / \vs\''il>dVg,= !^ e-''^'i^{-i{t))dt. 

We remark that a similar Gaussian beam quasimode construction was 
used to deal with partial data inverse problems in the paper [20] which 
was in preparation simultaneously with this manuscript. It is an interesting 
question whether other quasimode constructions could be used to extract 
more information about the potentials via Theorem 11.31 In particular, the 
following question is of interest. (By Theorem ll.4l we know that this question 
has a positive answer if A = for any (Mo,5'o) in which the ray transform 
is injective; on the other hand having A 7^ might help.) 

Question. Let (Afo, 50) be a compact oriented manifold with smooth bound- 
ary, and let A G M. Under which conditions on (Mq, 50) is the set M\ dense 
in the set of all bounded measures on Mq? 

The previous results are all based on extensions of the complex geomet- 
rical optics method. In the final results of this paper, we will use a com- 
pletely different approach and reduce the anisotropic Calderon problem to 
an inverse problem for the wave equation. To motivate this, note that the 
Laplace-Beltrami operator in a product type manifold (M x Mq, 5), where 
g = e® go and we now write t for the Euclidean variable, has the form 

By formally complexifying the t variable by 1 1— )• it (Wick rotation), we arrive 
at the wave operator 

Let us next describe a standard inverse problem for the wave equation. 
If (Mo,fifo) is a compact oriented manifold with smooth boundary, if G 
C(Mo), and if T > 0, consider the initial-boundary value problem 

(^2 - Ag„ + qo)u = in (0, T) x Mq, 

n(0) = dtu{<d) = 0, 

'"l(0,T)x9Mo = /• 

This problem has a unique solution u S C°°((0,r) x Mq) for any / G 
C^((0,T) X OMq), and we can define the hyperbolic DN map 

AfoX • ^r((0,T) X dMo) ^ C°°((0,r) X dMo), f ^ dM(o,T)xaMo- 
The inverse problem is to determine the metric go up to isometry and the 
potential qo from the knowledge of the DN map A^^gQ. This problem is 
closely related (and often equivalent) to an inverse boundary spectral prob- 
lem [18], to a multidimensional Borg-Levinson theorem [31], and also to an 
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inverse problem posed by Gel'fand [llj. In this paper, the wave equation 
inverse problem will be called the Gel'fand problem. 

The Gel'fand problem in the above formulation has a positive answer, 
under the natural necessary condition that T > diam(Mo,5o)- This follows 
from the boundary control method introduced by Belishev [3j and later 
developed by several authors; we refer to the book [I8j for further details. 
The boundary control method is based on three components: 

(1) Integration by parts {Blagovestchenskii identity): recover inner prod- 
ucts of solutions at a fixed time from the hyperbolic DN map. 

(2) Approximate controllability based on the unique continuation theo- 
rem of Tataru [l3]: solutions u{tQ, ■ ) are Lp' dense in the appropriate 
domain of influence. 

(3) Recovering the coefficients: this uses a boundary distance represen- 
tation of (Afo,5'o) together with projectors to domains of influence 
and special solutions such as Gaussian beams. 

An elliptic analogue of the Gel'fand problem is given by the following 
version of the anisotropic Calderon problem. Let (Mo,5'o) be a compact 
oriented manifold with smooth boundary, let go G C°°(Mo), and let T = 
M X Mq be an infinite cylinder equipped with the metric g = e (B go- Write 
{t,x) for the coordinates in M x Mq. Let also Spec(— A^q -|- qq) = {Xj}'jl^ 
where Ai < A2 < . . . is the set of Dirichlet eigenvalues of — A^,, -|- qq in 
(Mo, 50)- Consider the Schrodinger equation in T, 

{-dt - + Qq - \)u = in T, u\dT = /• 

For simplicity, in this introduction we assume that A G C \ [Ai,oo), that 
is, A is outside the continuous spectrum of — A^ -|- go in T (see Section [6] for 
the more general case A S C \ Spec(— A^^ -|- go))- Then for any / G C^{dT) 
the above equation has a unique solution C°°{T) n H^{T), and there is a 
linear DN map 

The next result shows that one can reconstruct the isometry class of an 
unknown manifold (Mo, 50) and also a potential go from the knowledge of 
dMo and the DN map A™gy(A). 

Theorem 1.5. Given the data {dT,A™g^{\)) for a fixed A G C \ [Ai,oo), 
where dT = M x dMo and A™g^(A) : C^{dT) C°°{dT) corresponds to 
the Schrodinger operator — A^-l-go on T, one can reconstruct a Riemannian 
manifold {Mo, go) isometric to (Mo,g'o) and the potential go. 

We obtain a uniqueness result as a consequence (Ai is the first Dirichlet 
eigenvalue of — A^p -|- go in Mo): 

Theorem 1.6. Let (Mo, 50) and (Mo, 50) be two compact manifolds with 
boundary ^Mo, and let go,go G C°°(Mo). If 

Af„"^„(A) = Af„"^-„(A) for some A G C \ ([Ai, 00) U [Ai, oo)), 

then go = ipQgo for some diffeomorphism ^po : Mo — )■ Mo with V'olaMo = Id, 
and also go = go- 
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As mentioned above, the proof involves a reduction from the eUiptic DN 
map to the hyperbohc DN map and the boundary control method. We also 
use the elliptic DN map on the transversal manifold, defined for A outside 
Spec(-Ago + qo) by 

^mlqoi^) ■ v\dMo ^ duv\9Mo, (-Ago + qo- X)v = in Mq. 
The argument proceeds roughly as follows: 

(1) Extend A^^JA) to act on weighted Sobolev spaces on dT. 

(2) If /c> 0, obtain A^^^^iX - k'^) for any h E C~(aMo) via 

(3) Recover Ajj"g^(/x) for /i G C from {Ajj^gjj(A— A;^)}fc>o by meromorphic 
continuation. 

(4) Recover A^^g^, from {K^^^ ^^^{^)} by Laplace transform in time. 

(5) Use the boundary control method to determine (Mq, (70) up to isom- 
etry and go from A^y^Q. 

It was proved in [19J that knowing the transversal DN maps {K^^^^ (/Li)}^gc 
is equivalent to knowing the DN map for the following equations: 

• Wave equation (0| — A^^ + qQ)u = in (0, 00) x Mq, 

• Heat equation {dt — Ag^ + qQ)u = in (0, 00) x Mq, 

• Schrodinger equation {idt — A^q + go)"" = in (0, 00) x Mq. 

Our results show that the elliptic equation [—df — A^^ + qQ)u = in M x Mq 
can be added to this list. 

Note that Theorems 11.51 and 11.61 are valid for arbitrary transversal mani- 
folds Mq without any restriction on the geometry, and they allow to recover 
both the transversal metric and the potential from the elliptic DN map. 
They are also the first uniqueness results for the Calderon problem that we 
are aware of which employ control theory methods (in particular approxi- 
mate controllability based on unique continuation for the wave equation). 
At the moment we can only show these results by going through the wave 
equation. It would be interesting to understand if there is a proof that would 
work with the elliptic equation directly. 

However, there is a severe restriction: the potential go has to be indepen- 
dent of the t variable, unlike in Theorems 1 1 . 1I[L3] where the scalar coefficient 
may depend on the Euclidean variable. In fact, the analogue of Theorem 1 1.6 1 
on a fixed compact manifold (M, g) CC (M x Mq, e©g'o) with two potentials 
independent of the t variable can easily be reduced to standard boundary 
determination results [H Section 8]. Of course, in the infinite cylinder T 
boundary determination is not so helpful and we use a reduction to the 
wave equation instead. 

The Wick rotation t ^ it suggests that the potential go should indeed be 
independent of the t variable, or at least real analytic in t, for this reduc- 
tion to the wave equation to work. The boundary control method for the 
wave equation also requires the coefficients to be independent of the time 
variable, although a variant of this method due to Eskin [10] allows lower 
order coefficients that are real analytic in time. The Gel'fand problem for 
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time-dependent coefficients is interesting in its own right; see [33], [35], [38] 
for some results when the background metric is Euchdean. 

This paper is structured as follows. Section [J is the introduction, and Sec- 
tion [2] gives the construction of complex geometrical optics solutions based 
on quasimodes and proves Theorem 11.31 Section [3] contains a direct con- 
struction of Gaussian beam quasimodes and the proofs of Theorems 1 1 . H [T7^ 
and II. 4[ In Section H] we give an alternative construction of Gaussian beam 
quasimodes based on a microlocal reduction via Fourier integral operators. 
The Calderon problem in an infinite cylinder is considered in the last two 
sections. Section [5] discusses the case where the spectral parameter is out- 
side the continuous spectrum and gives the proofs of Theorems 11.51 and 11.61 
and Section [6] extends these results to the case where the spectral parameter 
may be in the continuous spectrum but not in the set of thresholds. 
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M.S. are supported in part by the Academy of Finland (Finnish Centre 
of Excellence in Inverse Problems Research), and M.S. is also supported by 
an ERG Starting Grant. The authors would like to thank the organizers 
of the MSRI program on Inverse Problems and Applications in 2010, the 
Isaac Newton Institute program on Inverse Problems in 2011, the Fields In- 
stitute program on Geometry in Inverse Problems in 2012, and the Institut 
Mittag-Leffler program on Inverse Problems in 2013 where part of this work 
was carried out. D.DSF. would like to acknowledge the hospitality of the 
University of Jyvaskyla. 

2. Complex geometrical optics 

In this section we explain the construction of complex geometrical optics 
solutions based on quasimodes in (Mq, go) and use this construction to prove 
Theorem 1 1.3[ The argument is close to [SI Section 5]. 

We win assume that {M,g) is CTA with {M,g) CC (M x Mo,g) where 
g = c(e go), and {Mo, go) is a compact {n — l)-dimensional manifold with 
boundary. Let also q G L°°(M). We first note the identity 

c 4 (-Ag+q){c i u) = {-Ag + q)u 

where 

n-2 _ n-2 

g = e®go, q = c[q - c * Ag[c * )). 
This shows that it is enough to construct solutions to (— -|- q)u = 0. 

Writing x = {xi,x') for coordinates in M x Mo, the function ip{x) = xi is 
a limiting Carleman weight in a neighborhood of M [8j. In particular, we 
have the following solvability result which follows from [HI Section 4] (see 
also |2H Section 4] where one obtains solutions). 

Proposition 2.1. Let q E L°^{M). There exists tq > 1 such that whenever 
I ''"I ^ To, then for any / G L'^{M) the equation 

e™n-^9 + 9)e~™'r- = / inM 
has a solution r G H^(M) satisfying the estimates 

\\r\\HHM)<C\Tr~^\\f\\L^M), 0<a<l. 
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Consider complex frequencies 

s = T + iX, T real with |t| large, A complex and fixed. 

We are interested in finding complex geometrical optics solutions to the 
equation (— Ag + q)u = in M , having the form 

u = e~^^^{v + r). 

Here v = Vg will be an amplitude type term, and r = is a correction term 
with 1 1 r 1 1 (j^,/) — )• as |r| — )• oo. Further, we require certain asymptotic 
properties of as |t| — )• oo when A is kept fixed. 

A function u of the above type is a solution provided that 

where the right hand side is 

The point is to choose v so that ||/||l2(jv/) does not grow when |r| — oo, and 
to choose r so that e~*^^ir is the solution given by Proposition 12. li 

At this point we use the product structure on (M x Mq, g) where g = e© (7o> 
which implies that = 5^ + A^^. Consequently 

e--i(_Ag + q)e-'''H = {-df + 2sdi - - \, + q)v. 

This expression simplifies if we choose v independent of xi, that is, v = v{x'), 
and in this case 

Now ||/||l2(jv2) will not be too large with respect to |t| if f = Vs{x') is 
a quasimode or an approximate eigenfunction in the transversal manifold 
(Mo, 50), in the sense that 

||(-Ago - s^)vs||l2(Mo) = o(|r|), \\vs\\l2(Mo) = 0{l) 

as |t| — )• 00. 

The following result describes the complex geometrical optics solutions. 

Proposition 2.2. Let q G L°°{M), let tq be sufficiently large, and let A be 
a fixed real number. Suppose that {vs ; s = r + iA, |r| > tq} is a family of 
functions in LP'{Mq) satisfying 

||(-Ag„ - s^)vs\\l^Mo) = od"^!)' II^^IIl2(Mo) = 0(1) 

as |t| — )■ 00. Then for any r with |r| > tq there is a solution u G H^{M) of 
(— Ag + q)u = in M having the form 

where ||?'s||l2(jv2) = o(l) as |r| — )• 00. 

Proof. We first produce a solution of the equation (— A^ + q)u = having 
the form u = e~^^'^(fs + r^) as in the preceding discussion, and then define 

n — 2 

u = c ~u to obtain a corresponding solution of (— A^ + q)u = 0. □ 
The next result is slightly more general than Theorem II. 3i 



CALDERON PROBLEM IN TRANSVERSALLY ANISOTROPIC GEOMETRIES 11 



Proposition 2.3. Let {M,g) CC (M x Mo,g) be a CTA manifold, where 
g = c(e © (70)) and let gi, (?2 £ C(M). Let Ai, A2 G C, let {Tj)J^^ be sequence 
of positive numbers with — )• 00, and let 

Sj = Tj + iXl, tj = Tj + 1X2- 

Suppose that {vsj), (wtj) C C^(Mo) are sequences satisfying 

\\i-^90 - s'j)Vs,\\L2(Mo) = O(rj-), \\Vs,\\mMo) = 0(1), 
IK-Ago - t^)u;4j|i2(A^Q) = 0{Tj), Wwt^h^Mo) = 0{1) 
as j — )• 00, and in the weak topology of measures on Mq, 

lim Vs.wTdVgg = fix^M 

for some bounded measure /^Ai,A2 on Mq. If 



then 



Mo 



00 

e-^(>^^+^^)-^(c{qi-q2)){xi,x')dxj 



—00 



Here qi — q2 is extended by zero to M x Mq . 

Proof. We use Proposition 12.21 to find solutions of (— + qijUg^ = and 
(-Ag +q2)ut^ = 0, of the form 

where H^'sj ||l2(m)) Iktj llL2(Af) = o(l) as j — )• 00. Note that ut^ solves the 
equation (— A + q2)utj = in M. 

Next follows the usual integration by parts: we have 

/ {qi - q2)Us,Utj dV = [{AgUs^)utj - Us.i^gUt^)] dV 

Jm Jm 

[{duUs^)utj - Us^{duUtj)] dS 



IdM 

where the normal derivatives of the solutions are interpreted in the weak 
sense as elements in H~^/'^{dM). Using the condition C^^g^ = Cg^q^, there 
is some u £ H^{M) with (— A^ + q2)u = in M and 

u\dM = UsjldM, dyU\dM = dyUsj\dM- 

This shows that 

/ {qi - q2)us,ut^ dV = [{di,u)utj - uidiyUt^)] dS 
Jm JdM 

[{AgU)Ut^ -u{AgUt^)] dV 

M 



/ (92 - q2)uut^ dV = 0. 
Jm 
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Substituting the forms of the solutions Usj and ut^ in the last identity, we 
see that 

{qi - q2)e-'^^'+^''^'''c-'^Vs^W-dV = o(l) 

M 

using the norm estimate for the correction terms r^. and rt and the L? 
estimates for Vs^ and wt^- We now extend qi — q2 by zero to M x Mq and 
note that dVg{x) = c^^'^ dxi dVgg{x'). Then, taking the limit as j — )• oo 
and using the assumption that Vs uh' converges in the weak topology of 
measures, we obtain 



J Mq L-'-OO 



dfix^,X2{x') = 0. 



To be precise, we would like that the expression in brackets is a continuous 
function with respect to x' in Mq in order to take the limit. However, the con- 
dition Cg,gi = Cg^g2 implies by boundary determination that qi\dM = Q2\dM, 
and thus the zero extension of qi — q2 is in fact a continuous compactly 
supported function in M x Mq. The boundary determination result is essen- 
tially contained in [1,3', Proposition A.l] for the case n = 2, and a similar 
argument works also for n > 3 (see [U Section 8] for the case of DN maps 
with smooth qi and 52)- D 

Proof of Theorem This follows from Proposition 12.31 by taking Ai = A2 
real and taking Vs^ = Wg^ ■ □ 

At this point it is useful to compare the solutions in Proposition 12.21 to 
the ones appearing in [8l Section 5], where the additional assumption that 
(Mo,(7o) is simple was imposed. The complex geometrical optics solutions 
in [S], satisfying (— + q)u = in M, have the form 

Here V' is a real function chosen as a solution of an eikonal equation, and 
the amplitude a solves a complex transport equation in M. Since (Mo,5o) 
is simple these equations can be solved globally in M, and in fact ip only 
depends on x'. Then e~^'^^a satisfies 

as r — 7- 00. If a would be independent of xi, then e~'^'^^^^'^a{x') would be 
an approximate eigenfunction in Mq in the sense that 

(-A,„-r2)(e— '^a) = 0^.(Mo)(l). 

However, such functions are not quite sufficient to prove uniqueness results 
for the inverse problem. In |8j one instead employed amplitudes of the form 
a(xi,x') = e~^^^^d{x') which allow to exploit the Fourier transform in xi. 

There are two differences between Proposition 12.21 and the construction 
in [8], although the two are very closely related. The first one is that we use 
large complex frequencies s = T+i\ instead of large real frequencies r, which 
amounts to incorporating the factor e~*^^i from the amplitude a as part the 
complex frequency (thus making it possible to use the Fourier transform 
in xi). The second difference is roughly that instead of using approximate 



CALDERON PROBLEM IN TRANSVERSALLY ANISOTROPIC GEOMETRIES 13 



eigenfunctions e"*'^'^^^ ^a{x') with real frequency, we consider more general 
approximate eigenfunctions Vs{x') with slightly complex frequency. This ap- 
proach loses some generality since Vs is not allowed to depend on xi, but has 
the benefit that one can use much more general approximate eigenfunctions 
Vs{x') than those of the form e~*'^^^^'-'a(x') obtained from a global WKB 
construction on Mq. 

We end this section by noting that in Section 5 we will show that the 
semiclassical limit measures arising in Theorem 11.31 are invariant under the 
geodesic flow in a suitable sense. To do so, we will lift the measures associ- 
ated with the quasimodes v-r+ix to the cotangent bundle T*Mq and introduce 
the corresponding semiclassical measures. 

3. Gaussian beam quasimodes 

We will now give the Gaussian beam construction of approximate eigen- 
functions, or quasimodes, with desirable concentration properties. In fact, 
these quasimodes will concentrate near a geodesic in the high frequency 
limit. On a compact manifold without boundary, it is well known that one 
can find quasimodes concentrating near a stable closed geodesic for large 
real frequencies. We refer to |451 Section 10] and the references therein. 

The setup here is more flexible since there are no boundary conditions or 
global conditions on a closed manifold required of the family {vg}- Therefore, 
a construction of local nature is sufficient. We will give a direct argument 
analogous to the construction of Gaussian beams, which are approximate 
solutions of the wave equation localized near a geodesic [18] . The fact that we 
need approximate eigenfunctions with slightly complex frequencies instead 
of real ones will not present any complications. A version of this construction 
that also takes into account possible reflections is given in |2Uj . 

For most of this section we will write {M,g) for the transversal mani- 
fold instead of (-/Vfo,(?o) in order to simplify notation. Let {M,g) be an m- 
dimensional compact oriented manifold with smooth boundary (thus m = 
n — 1 > 2). Recall that a unit speed geodesic 7 : [0, L] — M is called non- 
tangential if 7(0), 7(L) are nontangential vectors on dM and 7(i) G M™* 
for < t < L. Theorem 11.41 is the following statement. 

Proposition 3.1. Let 7 : [0, L] — t- M be a nontangential geodesic, and let 
A G M. For any K > there is a family of functions (vs) C C°°(M), where 
s = T + iX and r > 1, such that 

\\{-Ag - S^)Vs\\mM) = 0{t-^), \\Vs\\l2(Mo) = 0(1) 
as r —7- 00 and for any -0 G C{M) one has 

lim [ \vs\'^ijdVq = [ e-2^V(7(*)) rfi- 
^^'^ Jm Jo 

In the case where (M, g) is simple, the method in [8] (although it was not 
written exactly in this way) reduces to using approximate eigenfunctions 
of the above type to recover attenuated geodesic ray transforms of desired 
quantities. In fact, a version of Proposition 13 . 1 1 on simple manifolds follows 
easily from the methods in [8]. 
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Proposition 3.2. Let {M,g) be simple, let A be a fixed real number, and 
let 7 : [0, -L] — )■ M be a nontangential geodesic. For any < a < 1 there is 
a family {vg ; s = t + iX, r > 1} in C°°(M) such that 

||(-A, - s')vs\\mM) = 0(r"), ||^;.||l2(m) = 0(1) 
as r — oo, and for any ip £ C{M) 

\vr+i\\^iljdVg ^ I e~'^^^'ip{'y{t))dt as r oo. 
iM Jo 

Proof. One first embeds {M,g) in a slightly larger simple manifold {D,g) 
and considers polar normal coordinates (r, 9) centered at a point u G D\M. 
There exist u £ D \ M and ^ S™'"^ so that 7 is part of the geodesic 
r I— )• (r, ^o) in ^ (fmy nontangential geodesic in M arises in this way for 
some uj and 6q). By using a WKB ansatz and choosing suitable solutions of 
the eikonal and transport equations as in ^ Section 5], the quasimodes at 
frequency s = t + iX can be chosen as 

Vs{r,e) = eng{r,e)\-'/%{e) 

where br G C°°(S'™~^) is an approximation of the delta function so that 

II^t||l2{s™-i) = ll&T||iy2,oo(5m-i) = 0(r"), 

|6t-P dS weakly as r — 00. 

A direct computation shows the required norm bounds, and we have for any 

{v-j-^ixl'^tl^ dVg ^ / e~'^^^'ip{r,9()) dr as r — )• 00. 
M Jo 

□ 

We now move to the proof of Proposition 13.11 The main difference to the 
case where {M,g) is simple is that the quasimodes can not be constructed 
using the WKB ansatz by solving eikonal and transport equations globally 
in M. Instead, we follow the construction of Gaussian beams: the eikonal 
and transport equations are only solved to high order on the geodesic, and 
we employ a complex phase function with Gaussian decay away from the 
geodesic. The phase function will be obtained by solving a matrix Riccati 
equation. 

We first record a few elementary lemmas (for proofs see [201 Section 7]). 

Lemma 3.3. Let (M, g) be a closed manifold, and let 7 : (a, 6) — )■ M be a 
unit speed geodesic segment having no loops. Then there are only finitely 
many times t G (a, b) for which 7 intersects itself at 7(i). 

Lemma 3.4. Let be a C°° map from a neighborhood of (a, 6) x {0} in 
into a smooth manifold such that -F|(a,6)x{o} is injective and DF(t,Q) is 
invertible for t G {a,b). If [ao,6o] is a closed subinterval of (a, 6), then F is 
a C°° diffeomorphism in some neighborhood of [oq, ^q] x {0} in M"-. 

The next lemma gives a system of Fermi coordinates near a geodesic that 
will be useful for the construction of Gaussian beam quasimodes. If the 
geodesic self- intersects, one needs several coordinate neighborhoods. 
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Lemma 3.5. Let (M, g) be a compact manifold without boundary, and 
assume that 7 : (a, 6) — )• M is a unit speed geodesic segment with no loops. 
Given a closed subinterval [ao,6o] of (a, 6) such that 7|[ao^;,o] self-intersects 
only at times tj with ao < ti < . . . < tjy < bo (set to = and t^+i = bo), 
there is an open cover {iUj,(pj)}^J~J of ^{[ao,bo]) consisting of coordinate 
neighborhoods having the following properties: 

(1) (pj{Uj) = Ij X B where Ij are open intervals and B = B{0, 6') is an 
open ball in W^~^ where 6' can be taken arbitrarily small, 

{2) ipj{j{t)) = {t,0)fovtGlj^ 

(3) tj only belongs to Ij and Ij n 1^ = unless \ j — k\ < 1, 

(4) ipj = ipk on ipj^iilj n Ik) X B). 

Further, the metric in these coordinates satisfies g^^\^(^t) = '^"''" j dig^^\^{t) = 0. 

Proof. The proof is based on Fermi coordinates. Choose {v2, ■ ■ ■ ,Vm} in 
T^(^ao)^ such that {'y{ao),V2, ■ ■ ■ ,Vm} is an orthonormal basis of r^^^o)-^- 
Let Ea{t) be the parallel transport of Va along the geodesic 7. Since 7(t) 
is also parallel along 7, the set {^{t), E2{t), . . . ,Em{t)} is an orthonormal 
basis of r^(t)M for t G (a, b). 
Define the function 

F : (a, b) X M"-i ^ M, F{t, y) = exp^(,) 

Here exp is the exponential map in (M, g) and a, /3 run from 2 to m. Then 
F{t, 0) = 7(t) and (with the ath coordinate vector) 

— F(t,se,)|^^o = E^{t), g-/{t,0) = i{t). 

Thus is a C°° map near (a, 6) x {0} such that DF{t, 0) is invertible for 
t £ {a,b). 

In the case where 7 does not self-intersect, -F|(a,6)x{o} is injective and 
Lemma 13.41 implies the existence of a single coordinate neighborhood of 
7([ao, 60]) so that (1) and (2) are satisfied (then (3) and (4) are void). In the 
general case, by Lemma [331 the geodesic segment 7|[aQ^f,o] only self-intersects 
at finitely many times tj with ao < ti < . . . < tjy < bo ■ For some sufficiently 
small 6, 7 is injective on the intervals (a, ti—6), (ti— 2(5, t2—S), . . . , {t!\/—26, b) 
and each interval intersects at most two of the others. Restricting the map F 
above to suitable neighborhoods corresponding to these intervals (or slightly 
smaller ones) and using Lemma [3.4l we obtain the required coordinate charts 
with (fj = F~'^\uy 

It remains to check the form of the metric in the coordinates x where 
= t and x'^ = with a = 2, . . . , m. Since the set {j{t), E2{t), . . . , Em{t)} 
is orthonormal, it follows that 

Thus also digjkl-y^t) = 0. We compute 

dagjk = da{dj,dk) = (Va„5j,9fc) + {dj^Vo^dk) 
where V is the Levi-Civita connection. Since 'Vg.dk = Vgi^dj, we have 

^9a^ll7(t) = "^didaljit) = ^j{t)Eait)\'y{t) = 
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because Ea is parallel along 7. Thus i9a5'ii|-y(j) = 0, and also 
dagwl^it) = Va^ (9/3)1 ^(t) = r^/3(7(i)) 

where F^-^ are the ChristofFel symbols. Considering the geodesic z{s) = 
exp^(^^^{sa°' Ea{t)) for some constants a", so that z{s) is given in the x coor- 
dinates by z{s) = {t, so?, . . . , sa^), the geodesic equation 

implies that for all a, /3 = 2, . . . ,m and / = 1, . . . , m we have 

Ti^m)a^a^ = 0. 
Since a" were arbitrary and F^^ = T^^, we obtain 

rL/3(7W) = 0. 

Thus daOipl^^t) = 0. Finally, 
= 0. 

We have proved that Sj^-^'^l^^j) = for all k. □ 

Proof of Proposition \3.1[ We begin by embedding (M, g) in some closed 
manifold {M,g), and extend 7 as a unit speed geodesic in M. Let e > 
be such that j{t) e M \ M for t E [-2e, 0) U (L, L + 2e] (here we use the 
fact that 7 is nontangential) . Our purpose is to construct a Gaussian beam 
quasimode near 7([— e, L + e]). 

Fix a point po = 7(^0) on 7([— e, L + e]) and let (t, y) be coordinates near 
Po, defined in a set C/ = {{t, y) ; \t — to\ < S, \y\ < 6'}, such that the geodesic 
near pQ is given by F = {(t, 0) ; |t — to| < S}, and 

g"'\r = 6^', ^^g'''\^ = 0. 

Here we write x = (t, y) where t = xi and y = {x2, ■ ■ ■ , Xm)- (Of course we 
will later use the coordinates in Lemma [3. 51 ) We will construct a quasimode 
Vg concentrated near F, having the form 

Vs = e''^a 

where s = t + iX, and Q and a are smooth complex functions near F with 
a supported in {\y\ < 6' /2}. 
We compute 

(-A - s^)v, = f 

where 

/ = e^"®(s2[((dG, dQ) - l)a] - is[2{de, da) + {AQ)a] - Aa). 

Here, the (7-inner product ( • , • ) has been extended as a complex bilinear 
form to complex valued tensors. We first choose G so that 

{dQ,de) = 1 to iVth order on F. (3.1) 
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In fact we look for B of the form B = X^jLo ®i where 

|a|=i 

We also write g^^ = Y1i=q + ''tv+i where 



By the properties of our coordinates, ffl^ = 5^^ and g'^-^ = 0. 

Choose Bo(t) = t and Bi(t,y) = 0. With the understanding that j, k run 
from 1 to m and a, /? run from 2 to m, we have 

g^'^djQdke -l = {l+gl' + .. .)(! + dtOi + . . .)(! + ^tBs + . . .) 
+ 2{gl" + . . .)(! + dtQ2 + . . ■)idy^Q2 + ■ ■ ■) 
+ (^"'^ + gf + . . .)(a,,e2 + dy^e, + .. .){dy^Q2 + 5,,B3 + ...)-! 
= [2atB2 + Vj,B2-Vj,B2 + 52'] 

+ [2atBp + 2Vj,B2-Vj,Bp + ^tQjdSk 

j,k<p 



p-2 



A^+ll 



1=2 j+k=p+i-l 1=0 j+fe=p+2-i 

k>2 'i<j,k<p 

We want to choose B2 so that the first term in brackets vanishes. Writing 
®2{t,y) = ^H{t)y • y where H{t) is a smooth complex symmetric matrix, it 
follows that H should satisfy the matrix Riccati equation 

Hit) + H{tf = F{t) 

where F(t) is the symmetric matrix such that g^{t,y) = —F{t)y-y. Choos- 
ing H(to) = Hq where Ho is some complex symmetric matrix with Im(i?o) 
positive definite, it follows that the Riccati equation has a unique smooth 
complex symmetric solution H{t) with lin(H{t)) positive definite [18, Lemma 
2.56]. This completes the construction of B2. From the lower order terms 
we can find B3, . . . , Bjv successively by solving linear first order ODEs on T 
with prescribed initial conditions at to- In this way, we obtain a smooth G 
satisfying 

The next step is the find a such that, up to a small error, 

s[2{de, da) + (Ae)a] - iAa = to Nth order on T. 
We look for a in the form 

a = T~^{ao + s~^a_i + . . . + a-N)x{y / 
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where x is a smooth function with x = 1 for \y\ < 1/4 and X = for 
\y\ > 1/2. Writing 77 = A 0, it is sufficient to determine Uj so that 

2{dQ, dao) + rjao = to iVth order on T, 
2{d&, da-i) + r]a-i — iAqq = to iVth order on F, 

2{d@, (ia_7v) + rjO—N — iAcL-(N-i) =0 fo A^th order on F. 

Consider oq = oqo + . . . + aoAf where aoj(t, y) is a polynomial of order j in 
y, and similarly let ?? = 7?o + • • • + Vn- The equation for ao becomes 

2(1 + + . . .)(1 + 9*02 + . . .)(5taoo + 9taoi + . . .) 

+ 4{gl" + ...)(! + 5^62 + . . .)(9j/.aoi + dy^ao2 + • • •) 
+ 2((5"^ + gf + . . .)(9y„e2 + dy^es + .. .)(5,«aoi + 5j,,ao2 + • • •) 

+ {m + m + ■■ ■)iaoo + aoi + . . .) 
= [25^000 + Voaoo] + [29tOoi + 2Vy62 • Vj^aoi + ??oooi + maoo] + 

Here 

rioit) = Ae(t,0) = dyjH^p{t)yp) = tr H{t). 
We want to choose aoo so that the first term in brackets vanishes, that is, 

dtaoo + ^(triJ(t))aoo = 0. 

This has the solution 

u\ ~h ftr.^'^ H(s) ds /, \ 

For later purposes we choose the constant as 

^detIm(ff(to)) 
Co = , (3.2) 



\Jkrn-ie ^y^'dy 



We obtain aoij • • • > QoAf successively by solving linear first order ODEs with 
prescribed initial conditions at tQ. The functions ai,. . . ,ai\f may be deter- 
mined in a similar way so that the required equations are satisfied to A?^th 
order on F. This completes the construction of a. 

We have constructed a function Vg = e'*®a in U where 

@it,y)=t+^H{t)yy + @, 



a{t,y) = T 4 {ao + s a_i + . . . + s a-N)x{y/^'), 

n\ —h ft tr H(s)ds 

ao(t, 0) = Coe ^ Jto ^ ' . 



Here 6 = 0(1 yp) and and each Uj are independent of r. Also, / = 

(—A — s'^)vs is of the form 

/ = e-®T^(s2/i2a + shi + ... + s-(^-i)/i_(jv-i) - zs-^A(a_Arx(y/5'))) 
where for each j one has hj = to iVth order on F. 
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To prove the norm estimates for Vs in U, note that 

IgisOj ^ g-AReOg-rlme ^ ^-At g- irlm(//(i))yy g-AO{|j/|2) g-rO{|j/|3) ^ 

Here lm.{H{t))y • y > c\y\'^ for {t,y) G [/ where c > depends on Hq and 5. 
This impHes that for t in a compact interval, after decreasing 6' if necessary, 
we have 



This shows that 



\vs{t,y)\<T'^e--^'-\y\\{y/6'). 



||(-A - .^)..||^.(^) < ||r^e-i-H^r2|yr+^ + r-^)||^.(^) = O(r^) 

as r — )• oo. The norm estimates for Vg in ?7 follow upon replacing N by 
2K + 3. 

For later purposes we record an additional estimate: if [/ fl dM ^ 0, the 
fact that the geodesic is nontangential allows to write dM locally in the 
{t,y) coordinates as {{t{y),y); \y\ < e} for some smooth function t = t{y). 
By choosing 6' small enough, we then have 

\\vs\\h(SM) = [ \vsitiy),y)\'dSiy) < [ r"^e-^-l^l'dy 



m — 1 



\y\<£ 

0(1) (3.3) 



as r ^ oo. 

We will now construct the quasimode Vg in M by gluing together quasi- 
modes defined on small pieces. Let 7([— e,L + e]) be covered by open sets 
as in Lemma [33] corresponding to intervals I^^^ (with the same 
6' for each U^^^) such that one can find quasimodes in each U^^\ We first 
find a function vi^'^ = e*'^®*°'a('^^ in U^'^^ as above, with some fixed initial 
conditions at t = — e for the ODEs determining and a^^\ Choose some 
t'o with 7(to) G U^^^ n U''^\ and construct a quasimode vi^^ = e^'^^'^'a^^-' in 
[/(I) by choosing the initial conditions for the ODEs for Q^^^ and a^^^ at 
to be the corresponding values of 0^*^^ and a^^^ at tg. Continuing in this way 
we obtain vi^\ . . . , vi^\ Let {xjit)} be a partition of unity near [— e, L + e] 
corresponding to the intervals {I^-'^}, let Xj(t, y) = Xi(^) in C/^-'^ and define 

r 



j=0 



Note that the ODEs for the phase functions and amplitudes have the same 
initial data in U^^^ and in U^^~^^\ which shows that we actually have vi'''^ = 
vi'''^^^ in U^^^ n U^^~^^\ In particular, if pi, . . . are the distinct points 
where the geodesic self- intersects, if < ti < . . . < t^/ are the times when 
the geodesic self-intersects, and if Vi, . . . , Vr are small balls centered at pj, 
then choosing 6' small enough we have a covering 

supp(t;,) n M C (uf^iVj) U {u^^^Wk) 
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where, in each Vj, the quasimode is a finite sum 



7(i;)=Pj 

and in each Wk there is some l{k) so that the quasimode is given by 

This shows that bounds for Vg and (—A — s'^)vs in M fohow from the 
corresponding bounds for each vi^^ . 
We stih need to verify the hmit 

\vr+i\\'^ip dVg -?■ [ e~'^^^il){'^{t))dt as r cx) 
iM Jo 

for any tp £ C{M). By a partition of unity, it is enough to consider functions 

tp with ij} £ CciVjCiM) and tp € CdWkCiM) (thus may be nonzero on dM). 

Let us begin with the case where ip G Cc{Wk H M) for some k. Then Vs = 

e"®a where G = t + \H{t)yy + 0{\y\^) and a = {ao + 0{T-^))x{y / 5'). 

Let p = \g\^^'^- We have 

\vr+ix\'^i:dVg 

M 

/ e-im{^fW)yye^-i/'0(|s/|3)gr-io(|y|2)^ 

(|ao(t, r-i/2y)|2 + o(^-i))^(y/^i/25/)2^(^^ r-^/2y)/9(t, r-^/^y) dt dy. 

Since Im(i?(i)) is positive definite and 5' is sufficiently small, the term 
g-im(_f/(t))j/-y (dominates the other exponentials and one obtains 

lim / \vr+i\\^iljdVg 



¥00 



'M 

e-^'^^^^''^^yy dy) \ao{t,0)\'^i;{t,0)p{t,0)dt. 
Evaluating the integral over y and using that p{t, 0) = 1 gives 



L 







T— >00 



lim / \vr+i\\'^ipdVr 



a 

M 



e-^y^'dy] I e-''^' ^{t,0)dt. 

J Jo Vdet Im(i/(t)) ' 

Here ao{t,0) = ao{tQ,0)e ^ Jto . Now we use the fact in [181 Lemma 

2.58] that solutions of the matrix Riccati equation have the property 

detlm(i?(i)) = detIm(i^(^o))e-'4*'•^^(^W)''^ 



2Xt l«o(i,n^'' 



It follows that |ao(^) 0)P/\/clet Im(if(t)) is constant in time. The choice 
(j3.2p fixes this constant and proves the limit for ip S CciW^ H M). 
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Now assume that ip S Cc{Vj D M), so that 

7(ii)=Pi 

It follows that 

\vs\'= E bff + E 

7(ti)=7(t;)=pj 

The computation above gives the right limit for each term. Therefore, 

it is enough to show that limits for the cross terms vanish as r — )• oo. 

Since all self-intersections must be transversal, and since deW(7(t/))is the 
covector corresponding to ^{ti) with respect to the metric, we may assume 
(by decreasing the sets Vj in the original construction if necessary) that 
Re((ie(') - de('')) is nonvanishing in Vj if 7(t/) = -i{tv) = pj but / / The 
cross terms lead to terms of the form 

JVjHM JvjnM 
where (/> = Re(0(^) — 0^'')) has nonvanishing gradient in Vj, and w^'^'^ = 

gisIm(e('-))g-ARe(<I.{'-))^{r)_ ^^^^ ^^^^^ ^Yiat 

lim / e'^'^w^^'^^i^dV = 0, I ^ l\ (3.4) 
'^^^ JvjnM 

showing that the cross terms vanish in the limit. To show (j3.4|) . let e > 0, 
and decompose tp = ipi + ip2 where ipi G C'^iyj H M) (-^i may be nonzero 
on dM) and \\ip2\\L'^{yjr\M) < £• Then 

I / e^-*u;(')^?/'2 dV\ < 11^2 11^2 ||V'2 < e 
JvjnM 

since ||_l2 < ||l2 < 1. For the smooth part ipi, we employ a non- 
stationary phase argument and integrate by parts using that 

e'^^ = —Lie'^'''), Lw = (\d4>r'^d(t),dw). 

IT 

This gives 

JvjnM JdM ir\d(j)\^ 

+ — [ e'^^L\w^^^^7Pi)dV 

Since Ib^'^^ ||l2(9m) = ^(1) by (|3.3p . the boundary term can be made arbi- 
trarily small as r — 7- oo. As for the last term, the worst behavior is when the 
transpose L* acts on g«*im{©'''')^ and these terms have bounds of the form 

|||d(Im(e«))|7;W||i2||t;('')||i2||V'i||Loo. 
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Here |d(Im(©(')))| < \y\ if {t,y) are coordinates along the geodesic segment 
corresponding to v^^\ and the computation above for ||v^'^||i,2 shows that 

|||d(Im(e«))|t;«||,.2||t;('')||i2||Vi||L=o < r-'/^. 

This finishes the proof of ()3.4p . □ 



In the end of this section, we switch back to writing (Mo,(7o) for the 
transversal manifold. Instead of using injectivity for the attenuated ray 
transform (see [HI Section 7] and [36] for injectivity results), we will reduce 
matters to the unattenuated ray transform. 

Proof of Theorem \1.4\ This is exactly Proposition 13.11 □ 



Proof of Theorem Assume the conditions in Theorem 11.21 and write 
q = c{qi — (?2)- As discussed in the end of the proof of Theorem 11.31 we can 
extend q by zero to M x Mq so that the extension, also denoted by g, is in 
Cc(K X Mq). Now, the combination of Theorems 11.31 and 11.41 implies that 

K2A,7(t))e-2^*(it = (3.5) 

for any A G M and for any nontangential geodesic 7 in Mq, where 

q(2X,x')= / e-^'^'''q{xi,x')dxi. 



00 



Thus the attenuated geodesic ray transform of q{2X, •), with constant at- 
tenuation — 2A, vanishes over all nontangential geodesies in Mq. 

Assume now that the unattenuated ray transform in Mq (the case A = 0) 
is injective. Evaluating ()3.5p at A = shows that 

for all nontangential geodesies 7. Injectivity of the ray transform then gives 
that q{0, • ) = in Mq. Next we differentiate ()3.5p with respect to A and 
evaluate at A = 0, to obtain 



2||(0,7(t))-2tg(0,7(t)) 



dt = 0. 



But since q{0, • ) = 0, this implies the vanishing of the ray transform of 
^(0, •) and hence also the vanishing of ^(0, •) in Mq. Taking higher 
derivatives with respect to A in (j3.5p and continuing this argument implies 
that 

d \ ^ 

— ) g(0, x') = for ah x' £ Mq 

OA J 

and for all /c > 0. Using that q{- ,x') is analytic as the Fourier transform of 
a compactly supported function, we see that q{^i.,x') = for all ^1 G M and 
x' G Mq. Thus 5 = 0, or gi = ^2 as required. □ 

Proof of Theorem\r^ First note that Cg = {{f,Agf); f G H^/^{dM)], 
where A^ is the DN map 

Ag : u\qm ^ duu\QM, ^gU = in M. 
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If {M,gi) and {M,g2) are two CTA manifolds in the same conformal class 
with Cg-^ = Cg2, we write 92 = 9 and gi = eg where c is some positive 
function. Then 

Acg = Ag- 

Boundary determination [HI Proposition 8.1] implies that c\dM = 1 and 
duc\dM = 0. Using fH', Proposition 8.2], this implies the following identity 
for DN maps of Schrodinger equations in (M,g), 

where Qc = Acg{c ^). Since {M,g) is a CTA manifold and the ray 
transform in the transversal manifold was assumed to be injective, we can 
now use Theorem 11.21 to conclude uniqueness of the potentials, —cQc = 0. 
But this implies that 

n-2 ^ n-2 , 

Acg{c- — ) = in M, c-—\9M = 1. 

Uniqueness in the Dirichlet problem implies that c = 1 in M, which shows 
that gi= g2- □ 

In the next section, we outline an alternative method for constructing 
quasimodes concentrating near a geodesic. The method is based on mi- 
crolocal reductions instead of the direct construction that was given above. 



4. MiCROLOCAL CONSTRUCTION 

Another possible approach to constructing quasimodes is a microlocal 
one; canonical quantization by a Fourier integral operator allows one to 
reduce the semiclassical operator Ag + to a simple form and construct 
the corresponding quasimodes. It will be convenient to use semiclassical 
conventions, and choose h = small parameter. We refer to [30j and 

|47j for a general presentation in semiclassical analysis. Let us nevertheless 
begin, for the convenience of the reader and to set our notations, by recalling 
a few definitions and results, which we will need in our exposition. 

4.1. Elements of semiclassical analysis. Semiclassical Sobolev spaces 
-ff^j.[ on a closed Riemannian manifold (or in Euclidean space) are defined like 
classical Sobolev spaces but are endowed with the following norms depending 
on the semiclassical parameter /i G (0, 1], 

\\n\\Hi^^ = (j:\myu\\iy 

where V-' are covariant derivatives on the Riemannian manifold. Semiclas- 
sical symbols of order k on T*W^ are smooth functions a on M^"^ depending 
on a parameter h G (0, 1] for which for all multiindices (a, /?) S N x N we 
have 

C«/3= sup (l + |C|')^|5°afa(x,C,/i)| <oo. 

(x,OeT*R™-,/ie{o,i] 
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The linear space of such symbols is denoted S^^y{T*W^). Pseudodifferential 
operators are defined through the semiclassical Weyl quantization 

(OVhaHx) = {^^h)-^JJ ei(--y>^a(^^,^,hy{y)dydC 

of a symbol a G 5^^j(r*M™') and we denote ^^^.[(M'") the corresponding space 
of operators. Symbols on the cotangent bundle on a compact manifold are 
smooth functions on T*M x (0, 1] which after cutoff to a coordinate patch 
pull back under local coordinates to symbols on T*W^. Pseudodifferential 
operators of order k on & compact manifold M are operators Ah : C°°{M) — )■ 
C°°{M) such that for all pairs of coordinate patches U,V and all cutoff 
functions if £ C^iU),^ G C^iV) 

• \\ipAhip\\j^,j^-N jjN) = 0{h°°) for all integers N if the supports of 

^ scl ' scl' 

and are disjoint, 

• ipAfiip written in local coordinates is a pseudodifferential operator 
Op;, a on with symbol a G S'^(T*M'"). 

We write ^'gpj(M) for the linear space of semiclassical pseudodifferential 
operators of order k on M. Using a partition of unity and local coordinates, 
it is possible to quantize any semiclassical symbol a G S^^^{T*M) into a 
pseudodifferential operator Op/^a G '^^^■^{M). Conversely, one can define a 
map which to any pseudodifferential operator A^ G ^^^.^(M) associates a 
class [a] of symbols in S^^^{M) / hS^^^ (M) called the semiclassical principal 
symbol of A such that Ah — Op^ a G h^^~^^{M). As usual one identifies a 
class of symbols with any of its representatives. 

Definition. A family u = {uh}o<h<ho of distributions on a closed compact 
manifold M or on is said to be tempered if there exists an integer N 

such that llw/tllji^-iv = 0{h~^). 

scl 

The semiclassical wavefront set WFsci(tt) of a tempered family of distri- 
butions u = {uh\o<h<ho on a compact manifold M (resp. W^) is the com- 
plement of the set of points (xo,Co) € T*M (resp. T*W^) for which there 
exists a symbol a G S'^.j such that, for some constant c > independent of 
h, one has |a(xo,^o)| > c and 

||(0p;,aK||i2 =0(^-). 

If Ah = Op/j a, one traditionally denotes by WFsci(^/i) the essential sup- 
port of a, i.e. the complement of points (x,^) in the cotangent bundle for 
which d"d^a = 0{h°°) near (x,^) for all a, 13. 

In the Euclidean space M™, there is an equivalent definition involving the 
semiclassical Fourier transform 

Definition. A point (xo,Co) € T*M"^ does not belong to the semiclassical 
wavefront set of a tempered family of distributions on if there exist 
smooth compactly supported functions XjV' which equal 1 near xq, respec- 
tively ^0 such that 

^Thixn) = 0{h^). 
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Remark. From the previous definition, the behaviour of the semiclassical 
wavefront set with respect to tensor products is clear: 

WF,,i{u(S)w) = {{x,y,tr]) : {x,0 £ WFsci(n), {y,r]) G WFsdW}. 

Example 4.1. The following examples of semiclassical wavefront sets of 
functions in the Euclidean space will be useful to our purposes: both are 
easily deducible from the definition involving the semiclassical Fourier trans- 
form. 

1. Coherent states or wave packets 

WFsei ((7r/i)-Te-2kl^-^oP+t(x-xo)-€o) = {{xo,(o)} 

(This is example (i) page 195 section 8.4.2 in |47]). 

2. Smooth functions independent of the semiclassical parameter h : 

WFsci(ti) = suppu X {0} 
(see Remark (ii) page 195 section 8.4.2 in [17]). 

We also recall the action of semiclassical Fourier integral operators whose 
canonical relation is the graph of a canonical transformation. Fourier in- 
tegral operators are operators whose Schwartz kernels are semiclassical La- 
grangian distributions associated with a Lagrangian manifold. We will con- 
sider Fourier integral operators associated with a Lagrangian manifold which 
is the graph of a canonical transformation. We denote vri : T*M™' — t- the 
first projection. A Fourier integral operator of order k associated with the 
graph 

G = {(x,C,?(x,e)) :(x,OgV^} 
of a canonical transformation q : V ^ W between two open sets V, W of 
r*M"* is an operator which maps distributions on X = 'tti{V) to distributions 
on y = TTi (W) whose kernel can be written modulo a smooth function which 
is 0{h°°) as the sum of terms of the form 

Uh{x,y) = (2^/i)-™y" et(^(^'«)-^-«)a(x,e,/i)^ie 

where a £ S^{T*W^) and (f is a generating function of the canonical trans- 
formation We recall that a function ip : V ^ W is a generating function 
of G if 

G = {{x,dMx,0,di:^ix,C),0 ■■ {^,0 e v^} 

in other words, the relation between the canonical transformation ? and the 
generating function (p is given by 

a{x,dxip) = {di:ip,C)- 

For notational purposes, one needs to introduce the twisted relation 

G' = {{x,y,^, -rj) : {x,^,y,r]) € G}. 

Indeed, the semiclassical wavefront set of the kernel Uh is contained in G' C 
T*{X X Y). 

One denotes I^di^ ^ ^) G') the space of such Fourier integral operators. 
The reason for adopting this notation is that one abuses notations by iden- 
tifying Fourier integral operators with their kernels which are semiclassical 
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Lagrangian distributions on X x Y with semiclassical wave front set con- 
tained in the Lagrangian submanifold G' of T*{X xY). If one or the other 
of the sets X,Y is a manifold without boundary, then the previous form has 
to be understood in local coordinates in x or y. 

Lemma 4.1. Let Uh € -^i5:i("^ ^ ^^G') be a Fourier integral operator as- 
sociated with the graph G of a canonical transformation : V W 
(mapping distributions on an open set X = tti{V) C to an open set 
Y = TTi{W) C M), then the semiclassical wavefront set is transformed un- 
der the action oi Uh in the following way 

WF,,i{UhUh) C ?(WF,,i(n;,) n V). 

We will also need a semiclassical version of Egorov's theorem. 

Theorem 4.2. Let Uh € I^^^i^"" x M,G'), Vh G /scf(M x M™,(G-i)') be 
two semiclassical Fourier integral operators respectively associated with the 
graph G of the canonical transformation c; and the graph G~^ oic;~^, and A G 
^gpj(M) a pseudodifferential operator then VhAUh is a pseudodifferential 
operator in M with principal symbol x(?*o) where x is the principal symbol 
of the pseudodifferential operator VhUh G ^'^^.[(M). 

In the classical setting, this is Theorem 25.3.5 in [16], for semiclassical 
versions one can refer to [27\ Theorem 4.7.8] and [47\ Theorem 11.5]. 

4.2. Semiclassical defect measures. It is time to introduce the notion of 
semiclassical defect measures which our introduction evoked and which lifts 
the measure used in our proofs to the cotangent bundle. Let {M,g) be a 
compact Riemannian manifold with boundary. We refer to [6] for a survey 
on semiclassical measures and to [471 Section 5.3]. 

Definition. Let {vj)'jLi be a bounded sequence of functions on M and 
{hj)JLi a sequence of reals in (0, 1] (called a sequence of scales) converging 
to 0. There exist subsequences (vjj.)fcli5 {^jk)T=i ^^'^ ^ positive Radon 
measure fi on r*M'°* such that for ah a G C^(r*M'"*) 

lim / Ah Vj^vj;dV= / adfi 

where Ah is a semiclassical pseudodifferential operator with principal sym- 
bol a and parameter hj^,. Such a measure is called a semiclassical defect 
measure associated to the sequences {vj)'^^ and {hj)JL^i. 

We are interested in the semiclassical defect measures associated with 
our family of quasimodes (fsj)°?,i for a sequence Sj = + i\ with {hj)'jL^ 
converging to 

\\i-h]Ag - {1 + iXhj)^)Vj\\L2^M) = Oihj), \\Vj\\L2^M) = 0{l) (4.1) 

as j — oo. 

Lemma 4.3. All semiclassical measures associated to the sequence of quasi- 
modes (j4.ip are supported in the cosphere bundle S*M^^^. 
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This is a consequence of Theorem 5.3 in [47j since the semiclassical princi- 
pal symbol of —h'^Ag — (1 + iXh)^ is — 1. The adaptation to the manifold 
case is straightforward. 

Lemma 4.4. All semiclassical measures associated to the sequence of quasi- 
modes (|4.1|) satisfy the following transport equation 

\Hp)i2 = 4A/i 

where Hp is the Hamiltonian vector field of the symbol p{x,^) = |^|^ = 

Proof. Let a G C~(r*M™*) be real- valued, and choose tp G C;?°(M™*) which 
equals one on the projection of suppa on M'°*. Since a is real valued, the 
pseudodifferential operator Ah^ = Op/^^. a is self-adjoint and we have 

1 2 

— {[Ah^,-h'jAg]i;vj,'iljVj)^,^^^ = — Im( - h]Ag{i;vj),Ah^{i;vj))^,^^j^ 



4X{Ah^Vj,Vj)^,^^^^ + o{l). (4.2) 



The principal symbol of the commutator ihj^[A,—h'jAg] is the Poisson 
bracket {a, \^\g} = —Hpa therefore the left-hand side term equals 

( OPh, iHpa)vj,Vj) ^2 + 0{hj) 

Passing to the limit in the equality ()4.2p . we finally get 

/ {Hpo) dji = AX I adji 
Jt'M Jt*m 

which proves the claim. □ 

Remark. If we were considering the semiclassical defect measure ft associ- 
ated with quasimodes on a closed manifold M then the transport equation 
would imply = e^'^^fi where (jjt denotes the cogeodesic flow on {M,g). 

4.3. Microlocal quasimodes. As in the previous section, to simplify no- 
tations we will write {M,g) for the transversal manifold instead of (Mq, (7o)- 
Thus, let (M, g) be an m-dimensional compact oriented manifold with smooth 
boundary and let 7 be a non-tangential geodesic. Once again, we embed 
{M,g) in some closed manifold {M,g), extend 7 as a unit speed geodesic in 
M and let e > be such that -f{t) G M \ M for t G [-2e, 0) U{L,L + 2e]. 
We recall that h = is our semiclassical parameter. After factorization 
of the operator 

Ag + {T + iXf = T^h^Ag + (1 + iXhf) 



= -T^ [yj-h'^Ag + 1 + iXh^ i^^-K^Ag -1-iXh 
it becomes clear that one has to seek Vs such that 

-h^A„ - 1 - iXh)vs = Oih^^^), s = h-^ + iX. 

' ^ HUM) ^ 
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In fact, we will construct an 0{h°°) quasimode. First, we will need the 
following proposition from [16j, which is a global version of the microlocal 
canonical reduction of a principal type operator. 

Proposition 4.5. Let / be a compact real interval and F : / — t- S*M be 
a non-closed cogeodesic curve, let Em = (0, . . . , 0, 1) € M"*, one can find 
a neighborhood V of the segment A = {((xi,0),em) : xi G 1} C M^™ 
and a smooth canonical transformation q : V ^ ?(^) from V to an open 
neighborhood <i{V) of F(/) such that 

^((xi,0),e^) = F(xi), c;*{^-l){x,0 = 6- 

This is (a non-homogeneous version of) Proposition 26.1.6 in [16j applied 
to the symbol a = yy^^^O^i^ ~ The proof follows Hormander's book 
and we provide it in the non-homogeneous case. Note that bicharacteristic 
curves of a are cogeodesic in the cosphere bundle. 

Proof. Assume that G /. Using local coordinates (y, rj) in a neighborhood 
Wr{o) of F(0) G 5*M, one can complete 



= a{y,r]) = ^/g(y~r]) - 1, (y, r?) G VFr(o) 

into a system (x,^) of symplectic coordinates near (0,5^) by the Darboux 
theorem. The map 

X ■■ V(^o,s^) Wr(o) 

is a canonical relation from a neighborhood V(^o,em) (Oi^m) in T*W^ to a 
neighborhood VFr(o) of F(0) in T*M such that 

X*a = ei, x(0,em) = F(0). 

Note that 

X((xi,0),em) = F(xi) 

for all xi close to since both functions satisfy the same system of ordinary 
differential equations 

OXi 

0), e^) = Ui, X((xi, 0), e^)} = Ha o x((xi, 0), £„) (4.3) 

ox I 

with same initial condition. In order to prove the lemma, it suffices to extend 
X to a neighborhood V of the segment A, by taking the maximal solution 
of the system of differential equations 

^ = H o 

dxi " ° ^ (4.4) 
^?(0,x',e) = x(0,x',O. 

Note that xi — )• <?(x,^) are bicharacteristic curves of a, hence cogeodesic 
curves in the cosphere bundle flowing through V. Therefore it follows from 
the assumption that F is not closed that the extension ^ is a diffeomorphism 
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if V is small enough. Furthermore, it is a canonical transformation: if 
{X, E) = ?~^(y, rj) then we have 

HaXi = 1, HaXj = for j = 2, . . . , m, HaEj = for j = 1, . . . , m, 
from which we deduce by the Jacobi identity 

Ha{Xj, Xk} = Ha{Ej,Ek} = Ha{Xj,Ek} = 

and therefore these Poisson brackets are constant along integral curves of 
Hp, in particular 

{Xj,Xk} = {xj,Xk} = 0, {Ej,Ek} = {Cj,^k} = 0, {Xj,Ek} = {xj,Ck} = Sjk- 
Since ? extends x we have 

?((xi,0),em) = x{{xi,0),em) = T{xi) 
and since is a solution of (j4.4p we also get 

a(?(x,e)) = a(?(0,x',O) = a(x(0,x',e')) = 6- 
This completes the construction of the canonical transformation <f. □ 

Our choice for F is the cogeodesic curve in S*M which projects on the ge- 
odesic 7 in M and we take / = [—e,L + e]. The next step is the quantization 
of such a canonical transformation. 

Proposition 4.6. Let F : / — )• S*M be a non-closed cogeodesic curve, and 
let ? be the canonical transformation introduced in Proposition 14. 5i For all 
A G M, there exist semiclassical Fourier integral operators Uh G /£.j(M'" x 
M, G'), Vh G I^^yiM X R"^, (G-i)') associated with the graphs G, resp. 
of the canonical transformation <f, resp. such that WF^diUhYh — 1), 

WFsci(V/i?7/j — 1) do not intersect F(/), resp. A, and such that 

Vh {^-h?^g - 1 - U/i) Uh = (hDi - iXh) + Rh + Ah 

where Rh € h°°%^^{M) and Ah e is such that WFsci(Ah)nA = 0. 

Proof. There exis10 semiclassical Fourier integral operators Uh € /^^.[(IR™ x 
M, G'), Vh G /o i(MxM™, (G-i)') such that WFsei(f/^°F^O-l), WF,,i{V^U^- 
1) do not intersect F(/), resp. A. By Egorov's theorem one has 

V^ [^-h^Ag - l) = V^U'hhDi + hRl 
where i^JJ G *sci(^)- I* 

remains to improve the remainder R^, and this can 
be done by further conjugation by two elliptic pseudodifferential operators 
C/i, V^ E ^U^n such that y.if/i - 1, - 1 G vI,--(M-) and 

[hD,, Ul] + G /i-M/--(M-). 

"'^One chooses f/h to be non-characteristic near A x r(J) and the construction of Vh is the 
standard construction of a parametrix (see Remark on bottom of page 27 after Definition 
25.3.4 in |16 for the classical case, Theorem 11.5 in [37] for the semiclassical case). 
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This can be done by choosing Uj^ = Op/jO where a ~ ^j>o^'^'^j ^he 
asymptotic sum of a sequence of symbols (aj)j>o satisfying the recursive 
equations 

-rdxittj +raj = -rj-i 

where r is the principal symbol of R^, r_i = 0, and rj-i is a principal 
symbol of the operator 



Rj = h-^-'^([hDi,Oph{ao+- ■ ■+h^aj)]+hRlOph{ao+- ■ ■+h^aj)j G ^ 
This sequence of equations can be explicitly solved and the solutions 

ao = exp ^ y r{yi,x') dyi^ 



scl 



are symbols of order j. 

^h^h "h. — <'h "h 



Taking Uh = U?Ul and Vh = V^VP we finally get 



Vh^-h^Ag - 1 - iXh) Uh = {hDi - iXh) + {VhUh - l){hDi - i\h) + R^ 

with Rh = V,l[hDi, Ul] + RlUh G h°°-^-^{W^). Because of the wavefront 
set properties of f/^ and F^, the remainder term = {V^Uj^ — l){hDi—iXh) 
has a semiclassical wavefront set which does not meet A. □ 

Having reduced the operator, it is now easy to construct quasimodes for 
the simple normal form hDi — iXh; in fact, we may as well choose a solution 
of the equation [di + X)v = 0, and take as our quasimode the function 

Vs = Uh{H{xi)e-^^Wh{x')), \\wh\\L\w-^) = 0{\) (4.5) 

where Wh is smooth and where i7 is a smooth cutoff function supported in 
[— 2, +oo) which equals 1 on [— l,+oo) so that H{xi)e~^^^Wh{x') is an 
function. That this could indeed be a possible quasimode is a consequence 
of the following relation 



UhVh [\J-h?Ag - 1 - iXh) Vs = Uh {hPi - iXh){H{xi)e-''''^Wh{x')\ 

=-i{hd:c^H)e-^^lWh 

+ UhRhiHe-^'^'whix')) + UhAhiHe-^'^'whix')) 
which leads to the estimate 

\{0phX){^-h^^9-l-iXh)vs\\^^ 

scl 

< \\{OphX)UhAh{He-^-^Wh{x'))\\Hi+0{h^) 

scl 

when X ^ C^{V) is a symbol which equals 1 close to r(/). It follows from 

WFsci{UhAh{He-^''wh{x')) n r(/) c (?(WF,ei(A)))' n r(i) = 

that we have 

WiOphXWhAhiHe-^-^Whix'mHi^^ = 0{hn. 
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Since x is localized in a neighbourhood of the cogeodesic F, we need 
an additional estimate away from r(7); in order to have such an estimate, 
we must impose on our quasimode that its semiclassical wave front set be 
contained in r(/). This means that we require 

WFsciK) = {(0,e„)}. (4.6) 

Lemma 4.7. The semiclassical wave front set of the quasimode Vg given 
by ()4.5p with the microlocal constraint (|4.6p is contained in the cogeodesic 
curve r(/): 

WFsciK-i+a)cr(7). 

Proof. Let V' bs a cutoff function. From the remark just after the defi- 
nition of wavefront sets involving the semiclassical Fourier transform and 
from the example 4.1.2, we deduce that the semiclassical wave front set of 
'4^{x)e~^^^Wh{x') is contained in the line A = {(j;i,0,em) £ IR^™' : xi G M}. 
The lemma follows then from Lemma 14.11 since the line A is mapped into 
the cogeodesic F by the canonical transformation <;". □ 



From the Lemma, we have WFsci ((y— /?Ag — 1 — iXh)vi^-i^ix) ^i^) 
and since 1 — x is supported away from F(/), we deduce 

'l-Op^x){^J-h^\-i-iXh)vs\\ =0{h^) 

scl 

Together with the previous estimate, this proves that Vs is a quasimode. 

Having constructed our quasimode, we proceed to the study of the corre- 
sponding semiclassical measure /x. Let a E C^(r*M'°*), we have 

/ (Op/, a)Vh-i+ix Vh-^+ix dV 
Jm 



U;,{Oph a)Uh{H{xi)e-^''^Wh{x')) H{xi)e-^^^Wh{x') dxi dx . (4.7) 

By Egorov theorem, the conjugated operator has a simple principal expres- 
sion 

Ul{Ovha)Uh = OpJx?*«) + hRh 

where x is the principal symbol of Uf^Uh and Rh G ^tc\- choose our 
function Wh to be a wave packet 

Wh = (7r/i)~^e~^l'''l^+^^'" 
and in the construction of one can take x to be one on F(/). 
Lemma 4.8. The semiclassical measure associated to the family 

V, = (7r/i)-'^/7(xi)e-^^^e-27^l^'l'+t^- 
is i?2(xi)e-2^^i dxi (g) 6^>=o^^=em- 
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Proof. Let a E C^iW"^), we have 
'l2(r™) 



We can take {x+y)/2 and (x — ?/)/2 as new coordinates and after integration, 
we get 



X e7r2^i«ie'^(l^'l'+l«'-""l')a(xi, x', ^i, C') da; dyi (iC- 
We let h tend to and obtain 

/oo 

which completes the proof. □ 

Using the lemma, Egorov's theorem and passing to the limit in (|4.7|) . we 
get 

rL+e 



T*M J-e 
L 



JO 



since x equals one on r(/). One can sum up our construction in the following 
theorem. 

Theorem 4.9. For any non-tangential geodesic on a compact Riemann- 
ian manifold (Mo,go) with boundary, there exists a family of quasimodes 

{vh-^+ix)he{o,i) such that 

ih'Ag, + {l + iXhf)vf,-^+,^ = 0{h^) 

lb/i-i+iAllL2(Mo) = 

with semiclassical wave front set contained in the cogeodesic T projecting 
on 7 and with associated semiclassical measure /i on Mq^^ given by 

adfx= [ a{T{xi))e'^''^ dxi. 



lT*Mo Jo 

From this alternate construction, one can also deduce Theorem 11.41 



5. Calderon problem in a cylinder 

In this section we will prove Theorems 11.51 and 11.61 which consider an 
inverse problem in the infinite cylinder T = M x Mq with metric g = e(B go- 
Here (Mo,go) is a compact oriented m-dimensional manifold with smooth 
boundary, m > 2. We write (t, x) for coordinates on T where t is the 
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Euclidean coordinate and x are coordinates on Mq. The Laplace-Beltrami 
operator in T is given by 

We consider more generally the Schrodinger operator on T, 

-A + qo = -5| - Ago + qo 

where go € C°°{Mo) is real valued. It will be crucial that the coefficients go 
and Qo are independent of the t variable. 

The first point is to set up boundary measurements related to the Dirichlet 
problem 

{-d^ - Ago + qo)u = in T, u = ^ on dT. 

The spectral properties of the Schrodinger operator in the infinite cylinder 
are different from those on a compact manifold because of the presence of 
continuous spectrum. Let Ai < A2 <•••—)■ cjo be the Dirichlet eigenvalues 
of -Ago + qo in Mq, write Spec(-Ago + go) = {Ai, A2, . . .}, and let {Mill be 
an orthonormal basis of L^(Mo) consisting of eigenfunctions which satisfy 
(-A3,, + qo)(l)i = Xi^i in Mo, G ^^oH^o)- 
We will see the following facts: 

• -A + qo with domain H'^{T) n H^{T) is self-adjoint on L^T), 

• the spectrum of —A + qo is [Ai, 00), 

• if A G C \ [Ai, 00) then for any 5 € M 

(-A + qo- A)-i : Lj{T) ^ {u e H]{T) ; u\aT = 0}, 

• if A e [Ai,cx)) and A ^ Spcc(— A^y + qo) then for any S > 1/2 the 
following limiting absorption principle holds: 

{-A + qo-X- iOr^ : Ll{T) ^ {u e HI^{T) ; u\aT = 0}. 

The case of thresholds (or resonances) , where A = A; , is special and will not 
be considered here. 

Let L2(r) = L2(T, dV) be the standard space in T, and let H'{T) be 
the corresponding L?' Sobolev spaces. Since Mq is compact, we define 

HLciT) = {f;fe H%[-R,R] X Mo) for any R > 0}. 

Writing (t) = (1 + 1^)^/^, we introduce for s > the weighted spaces 

Lj{T) = {feLUT); {t)'feL\T)}, 

HUT) = {fe HUT) ; {t)'f G H'(,T)}, 

Hlo{T) = {feHUT);f\&T = 0}. 

Also, H^{T) = {/ G H^{T) ; /|ar = 0}. We define, in the L'^{T) duality, 

H-\T) = {H^iT))*. 

If s > 1/2 define the abstract trace spaces 

H'{dT) = H'+^/^{T)/{H'+'^/^{T) n H^{T)), 

Hi{dT) = if|+^/^(r)/(iJ|+^/^(r) n h^{t)). 

Since dMo is smooth and compact, these spaces can also be identified with 
standard weighted Sobolev spaces on dT. 
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In this section we will assume that A is not in the spectrum [Ai, oo) (the 
general case A G C \ {Ai, A2, . . .} is considered in the next section). The 
following proposition shows that there is a well defined DN map ^'^Q^qoi^) 
related to the operator —A + go ~ ''^ in the cylinder T. 

Proposition 5.1. If A G C \ [Ai, 00), then for any / G H^I'^{dT) there is a 
unique solution u G (T) of the equation 

(-A + go-A)u = inT, u\ffr = f- 

If f e C^{dT) then u G C°°(r) and there is a hnear map 

Ajo,,o(^) ■■ CTidT) ^ c-(ar), / ^ dMdT. 

For any G M, this map extends as a bounded linear map 
Kl^^^{\):Hf\dT)^Hl'\dT). 

The first easy observation is that one has unique solvability of the Dirich- 
let problem in T for certain frequencies by the usual Lax-Milgram type 
argument. The key point is that dT is compact in the direction transverse 
to t, so it has a Poincare inequality in the standard Sobolev spaces. 

Lemma 5.2. Let A satisfy —00 < A < Ai. For any F G H~^{T) the 
equation (—A + qo — X)u = F in T has a unique solution u G Hq{T), and 
one has < C||F||^-i(2^). 

Proof. We consider the bilinear form 

B{u, = y {{dt,xU, dt,xv) + qouv - Xuv) dV, u,v e Hq{T). 
This is a bounded symmetric sesquilinear form on Hq{T) and satisfies 

/oo /■ 
/ {\dtu\'^ + \dxu\'^ + qo\u\'^ - A|u|2) dVgo dt. 
-00 J Mo 

We note the inequality 

/ {\dxv\^ + qo\v\^) dVg, > Ai / \v\^ dVg„ v G H^{Mo). 

This shows that B{u,u) > e\\du\\'j^2 (rp^ + (Ai(l — e) — £||(?o||l°° — -^)ll^llL2(y) 
for < £ < 1. If £ is sufficiently small we see that B is coercive and the 
Riesz representation theorem shows the existence of a unique solution of 
(-A + qo-X)u = F in T. □ 

Elliptic regularity for the previous problem, even in weighted spaces, can 
be proved by a Fourier analysis argument. 

Lemma 5.3. Let A G C \ [Ai, 00) and let 5 G M. For any F G Lj{T) there 
is a unique solution u G HIq{T) of the equation (—A + q^ — X)u = F inT. 
Further, u G H^{T) and ||«||ij2(r) < Q,A||i^||L2(T)- 
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Proof. Write v{t,l) = {v{t, ■)-,4>i)l'^{Mq) foi^ the partial Fourier coefficients. 
If F G L|(r) we note that 

-oo J AIq J —oo 

OO 

1=1 

This imphes that F{- ,1) G L'j{]R.) for all /, and the Fourier transform satisfies 
F{ - ,1) G H^{M.). Formally, the equation (—A + qo — A)tt = F reduces to a 
system of ODEs: 

{-df + Xi- X)ui-,l) = F{-,1) on R, for 1 = 1,2,.... 

By taking Fourier transforms in t (with dual variable rj), we obtain 

u{t,l) = |-^--l_^F(r,,/)| . (5.1) 

Uniqueness follows immediately since if n G {{^^{T) solves the equation 
(-A + qo- X)u = in r for some 5 G M, then u{- ,1) G LjiM) and by taking 
Fourier transforms (ry^ + A; — X)u{r],l) = for r/ G M and for all I. Here 
r/^ + A; — A is never zero using that A G C \ [Ai, oo), so u( • , /) =0 for all / 
and M = 0. 

We move to existence and let F G L'g{T). If u{t,l) is defined by ()5.ip . 
then for k > \5\ we have 

II^(-'0IIl2(r) = J)\\hS{r) < IK^^ + Ai- A)"^||^yfc,oo(]g)||F(-,Olliyi(R)- 

We need to estimate the W'^''^ norm uniformly in I, using the condition 
A G C \ [Ai,oo). Note that 

iri' + k- A)-i = -L (5.2) 

2zi \ri- zi ri + zij 

where = A — A^ = a — (A^ — Ai) + z/3 and a = Re(A) — Ai, /3 = Im(A) (we 
assume zi £ {z; Re(z) > 0, Im(z) 7^ 0} U {ir ; r > 0}). Write 

fit) = faAt) = - t + ^/?)'^')' t > 0- 

When a < we have |/(i)| > |/(0)| for t > 0, and when a > we have 
\f{t)\ > \f{a)\ fort>a and |/(t)| > Cc,,/^ for < t < a. This shows that 

|Im(z;)| > Ca,/3 > 0. 

Since 

and since = |A - A^] > |Im(A)| if Im(A) / and |A - A;] > |Re(A) - Ai| 
when Im(A) = 0, we have 

||(r/^ + A/ - A)~^||^yfc,oo(]R) < Ck,a,i3- 

Thus 

M-,l)\\Lj{R) < C'<5,a||^(-,0IIl2(K)- 
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We define 

TV 

UN{t,x) = y^^u{t,l)(pi{x). 
1=1 

It follows that when M <N, 

N N 
l=M+l l=M+l 

Thus (utv) is a Cauchy sequence in L'^iT) and converges to some u G L'j{T). 
Since 

N 

(-A + qo- X)uN{t, x)=Y, ^(*' ^)Mx), 

1=1 

we obtain that u is a distributional solution of (—A + qo — X)u = F in T. 

We next show that u G Hj{T). The expression for u{t,l) together with 
()5.2p implies that 

\\dtu{.MLi + \k-\\''^\\u{-MLi<c\\F{.Mq- 

Thus \\dtu\\L2 < C\\F\\l2 and 

o o 

|2 



l|Vgo^^iv(i)llL2(Mo) = i-^ao^N{t),UN{t))L2(Mo) 



N 



Yi^i - ^Mt, 01' + ((A - qo)uN{t),UN{t))mM,). 



1=1 



Finally, we check that u £ Hj{T). Note that 



Consequently ||Vg(,M||^2(j.) < C||F||2^2(y), and then also u € H^q{T). 
{Xi - X)uit, /) = -^-1 I ^j^F{r^, I) 



and so, for > 

\\{Xi - X)u{ ■ , 1)\\l2 < \\zf/{rf - )||;y.,o. ||F( • , OIIl|. 
We use ()5.3p to write 

d V _ CjZi {rj + ziy+^ -{rj- ziy+^ 



drij ri'^-zf 2 (r/2 - zfy+i 

We choose L = L{X,go) so large that for Re{zf ) < for I > L. Then 
\r]'^-zf\ > for / > L, and also \zf /{rj^ - zf)\ < Cx,go for / = 1, . . . , L - 1. 
Also |r//(r/2 - \zi/{r]'^ - zf)\ < Cx,go- It follows that 

\\{Xi-X)u{-,l)h2<C\\F{.,l)^. 
Therefore (— A^^ + qo — X)u G L'j{T). By elliptic regularity 

INlo^it, ■)\\lHMo) < C{\\u{t, ■)\\l2(Mo) + \Wou{t, ■)\\lhmo)) 
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for almost every t. Consequently V^^u E L'j{T), which hnplies that also 
dfu E L^(r). By a short argument we obtain V g^dtu E L'j{T) and the proof 
is finished. □ 

The previous lemma also implies self-adjointness. 

Lemma 5.4. If qq E C°°(Mo) is real valued, then the operator —A + qq 
with domain H'^{T) D H^{T) is self-adjoint on ^^(T). 

Proof. The operator —A + qq with this domain is densely defined and sym- 
metric, and by Lemma 15.31 the range of —A + qq ±i is all of L'^(T). □ 

Proposition [5TT] is a straightforward consequence of Lemma [5. 3 1 and elliptic 
regularity (see also the proof of Proposition 16.11 below). The key point in 
the proof of the uniqueness result. Theorem 11.51 is the following connection 
between the DN map for Schrodinger operator —df — A^^ + in T and 
the DN map for the transversal Schrodinger operator — A^^ + go in -^o- We 
define the transversal DN map at energy E C \ {Ai, A2, ■ ■ ■} as 

A*^°gJ/i) : H^/\dMo) ^ H'/\dMo), h ^ d,Vh\dMo 

where Vh is the unique solution of the Dirichlet problem 

{-^go+Qo- f^W = in Mo, Vh\dMo = h. 

Proposition 5.5. If A E C \ [Ai,oo) and if /c E M, then 

In particular, the expression on the right is independent of the t variable. 

Proof. Let h E H^/^{dMo), and let Vh E H^{Mo) solve 

{-A + qo- {X- k'^))vh = in Mo, Vh\dMo = h. 

Note that since A ^ [Ai, 00), the number X — is not a Dirichlet eigenvalue 
of —Ago + go and there is a unique solution Vh- Define 

f{t,x) = e'''h{x). 

Since k is real, we have / E H^^'^{dT) for any 6 < —1/2. The function 
u{t,x) = e^^^Vh{x) is in H^{T) and solves 

{-dt - Agg +qQ- X)u = in T, u\aT = /• 

Thus 

This proves the result. □ 

We can now prove Theorems 11.51 and 11.61 showing that the DN map 
Aj^ g^(A) at a fixed energy A ^ [Ai, 00) uniquely determines the metric go up 
to isometry and also the potential q^. 
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Proof of Theorem \1.5[ Suppose one is given the manifold dT = R x Mq and 
the map Aj^ ^^(A) : C^{dT) C°°{dT) for some fixed A € C \ [Ai,oo). We 

Q /O 

may assume that BMq is known. Since C^{dT) is dense in H^'^{dT) for all 
6, we also know the map 

Al^JX):Hf{dT)^Hy\dT) 

for all 5. 

bmce e*''*/i G Hf'^{dT) whenever G M and 5 < -1/2, we may compute 
the map 

^gtM - • H'^\dMo) ^ Hy\dMo) 
for all A; G M from the knowledge of {dT, A'^^ g^{X)) by Proposition 15.51 
Since ^ i— >• Ag^^^^ (/x) is a meromorphic operator valued function whose poles 
are contained in {Ai, A2, . . .} [18\ Lemma 4.5], this information determines 
A^°gg(//) for all jj, in the complex plane by analytic continuation. This is 
equivalent to knowing the DN map for the wave equation 9| — Ag^ + qq in 
Mq X {t > 0} [18\ Chapter 4]. The boundary control method then allows 
to construct a manifold isometric to {Mo,gQ) and the potential qo from the 
DN map for the wave equation. See [TBJ, [191 more details. □ 

6. CALDERON problem in a cylinder: CONTINUOUS SPECTRUM 

Assume the conditions in the first paragraph of Section [5l We will next 
consider the case when A is in the continuous spectrum [Ai,oo) but outside 
the set of thresholds {Ai, A2, . . .}. In this case the Schrodinger equation in 
T admits generalized eigenfunctions, and a radiation condition is required 
for uniqueness of solutions and for the definition of the DN map. 

Proposition 6.1. Let A G [Ai,cxd) \ {Ai,A2,...}, choose Zq ^ 1 so that 
Xig < A < Ai„+i, let 5 > 1/2, and let m > 2. Then for any Hj^'^^^idT), the 
equation 

(-A + go - A)u = in T, u\qt = f 
has a unique solution u G H'^^{T) satisfying the outgoing radiation condition 

{dt =F Xi)u{t, /) as t — > ±00 for all 1 < / < /q. 

If / G C^{dT), then u G C°°{T) and there is a hnear map 

^IaM) ■■ CTidT) ^ C^{dT), f ^ dMdT. 
For any 8 > 1/2, this map extends as a bounded linear map 
^l,,M)--K'"'\dT)^H"^~''\dT). 

Recall that when A G C \ [Ai,oo), the main point in the reduction from 
the Calderon problem in the cylinder to the boundary control method was 
Proposition 15.51 This result states that 

e»%^o,„(A - k^)h = Al^JX){e^'''h\sT), h G H^/\dMo). 

This identity does not directly generalize to the case where A is in the con- 
tinuous spectrum, because the boundary value e^^^h\gT on the right hand 

o /o 

side is not in {dT) for 6 > 1/2. However, by using suitable cutoff 



CALDERON PROBLEM IN TRANSVERSALLY ANISOTROPIC GEOMETRIES 39 



and averaging arguments we can still recover the transversal DN maps from 
Ajo,.o(A)- 

Proposition 6.2. Let A G [Ai,oo)\{Ai, A2, . . .}, let A: G M, and assmne that 
\-k^ i {Ai, A2, . . .}. There is a family {^r)r>x C Cfi^) with -^^{1) = 1 
for \t\ < R, such that 

e'^'^ZoiX - k')h = Jim l\l^^^{X){e'''^n'{t)h\aT) clR' 

pointwise on dT for any h £ C°°(5Mo). 

Applying this result instead of Proposition [^31 we obtain an extension of 
Theorem 11.51 to the case where A is in the continuous spectrum (but not in 
the set of thresholds) with exactly the same proof. 

Proposition 6.3. Given the data (ST, Aj^ ^^(A)) for a fixed A € [Ai,oo) \ 
{Ai,A2,...}, where dT = R x OMq and Ajj,,j„(A) : C^{dT) C°°{dT) 
corresponds to the Schrodinger operator —A + go on T, one can reconstruct 
a Riemannian manifold (Mo,5o) isometric to (Mo,go) and the potential go. 

We now move to the proofs of Propositions 16.11 and 16. 2[ The first step 
is an existence and uniqueness result for the inhomogeneous Schrodinger 
equation in the cylinder. 

Lemma 6.4. Let A E [Ai,oo) \ {Ai,A2,...} and choose Iq > 1 so that 
A^o < A < Xiq+i- Let (5 > 1/2, let /i G M, and let m > be an integer. For 
any F = Fi + F2 where Fi G H^{T), F2 £ HJ^{T) and 

lo 00 
Flit, x) = Y, Flit, F2it,x)= J2 F2it,l)Mx), 

1=1 l=lo+l 

there is a solution u = ui + U2 of the equation 

i-d^ - Ag, + qo - X)u = F inT 
where m E i/™+^(r) n Hl^QiT) and U2 G n i^^,o(^) are of the 

form m = Y!i=i uii- ,l)4>i and U2 = T.'Zio+i '"2( • ,0</'«- Further, 

lkl|l/f™+2(T) ^ C\\Fi\\h^(^t), \W2\\h]J'+^{T) ^ C'II-^2||h™{T)- 

The solution is unique up to an element of the form 

1=1 

where are constants. If one assumes the outgoing radiation condition 

idt =F i\/A - Xi)uit, /) as t — )■ ±00 for all 1 < / < Iq, 
then the solution u is unique. 

Proof. For uniqueness, suppose that u G H^QiT) for some real number r 
solves i—df — + go — X)u = in T. Then the partial Fourier coefficients 
■u(t, /) satisfy 

i-df + A)m( -,0=0 in R, for all I > 1. 
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If / > /q + 1) then A/ — A > and by taking Fourier transforms we see 
that the only tempered distribution solving this equation on M is zero. Thus 
u{- ,1) =0 for I > Iq + 1. If 1 < / < /q) then A/ — A < and an easy argument 
shows that the only distributional solution of the above equation for u{ - ,1) 
is 

u{t, I) = c+e'"^^^ + cle-'"^^^ 

for some constants cf. If the radiation condition holds it follows that = 0, 
concluding the proof of uniqueness. 

For existence, let first m = 0. We define 

n2(r?, = ^-^1 ^2 ^!^^_^ F2(^,0}, />/o + l. 

Since A/ — A > for I > + the proof of Lemma 15.31 shows that the 
function U2 = Yl'i^io+i^2{- J)4'i solves the equation {—d^ — Ag^ + qo — 
X)U2 = F2 and has the required properties. The function ui is obtained as 
ui = X^llx ■ ) where ui{- ,1) should satisfy 

{-df + \)ui{-,l) = Fi{-,1) in M, for 1 < / < Iq. 

We choose the solution 

= i?o(A-A, + iO)Fi(-,/) 

where Roiz) = {—d^ — z)~^ is the resolvent of the Laplacian on the real line, 
and Ro{s + iO) is the outgoing resolvent at energy s > 0. If > one has 
the well known formula (which follows from a direct computation) 

/oo • 
G{t-t')f{t')dt', G(i) = -e^^N. 

Agmon's limiting absorption principle (see [16, Section 14.3], [34, Section 
XIII. 8]) gives that ui{-,l) G if^^(M), and then ui is a solution of (— 9| — 
AgQ+qQ—X)ui = Fi with the required properties (also satisfying the outgoing 
radiation condition). This concludes the proof for the case m = 0. 

We show the case of general m by induction: the statement has been 
proved for m = 0, and we assume that it holds for values up to m — 1. Let 
F2 G i?™(T) (the proof for Fi is analogous). Then also F2 G HJ^-^T), and 
by the inductive hypothesis there is a unique solution u G HJ^^^ (T) of 

{-d'^ - Ago + go - A)u = F2 in T, = 0. 

Then u satisfies the same equation with right hand side F2 G H^~-' (T) , 
so we have d{u ^ Hj^^'^ ■' (T) for j > 1 with 

\\&lu\\^rn + 2-3^j,^ < C 1 1 ^2 1 1 //m (T) . 

The equation implies that for almost every t, 

-Ag^uit, • ) = {d^ - go + X)u{t, • ) + i^2(t, • ) in ^o- 
Since also u{t, ■ )\dMo = 0) elliptic regularity implies that 

\\u{t, • )||j/,n+2(Mo) < C\\{dt -qo + A)n(t, • ) + F2{t, • )IIh™(Mo)- 
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Multiplying by (t)'* and taking the norm over M, we obtain 

ll'"llL2(R;H'"+2(Afo)) < C'(II^IU™(T) + l|-^2|U™(T)) < C|| -^2 |U™ (T) • 

Thus u,dtu,V,j;U e HJ^+^{T), showing that u E HJ^+'^{T) with the right 
bounds. □ 

Proof of Proposition [6A[ Let 6 > 1/2, and suppose that / G -ff™ ^^^((?T). 
Choose Ef G i^f (T) with Ef\aT = f and ||^^|Uj.(t) < C||/||^^„-i/2(^^^. 
We look for a solution of 

{-A + qo-X)u = in r, u\qt = f 
having the form u = + w. Thus, we obtain the equivalent equation 

(—A + go — = F in r, w\qt = 
where F = —{—IS.+qQ-\)Ef. Since ||-F||^m-2j-y^ < C'||/||^m-i/2j.g^-|, Lemma 
3] shows that there is a unique solution w G H^g(T) n H\^q{T) satisfying 



(c?j =F \i)w{t, Z) — ;> as t — )- ±00 for ah 1 < ^ < /q. 

One also has (t) ^ C\\f\\ m-i/2,r,^.. Thus we have a unique solution 

It = E''^ + w ^ H^g(T) to the original problem, satisfying the same radiation 
condition as w since {dt =F — ( • i ^ -^5(1^) for all L We also have 



< C\\J \\^^-l/2^Q^y 

The result follows. □ 

Before the proof of Proposition 16.21 we record some further properties of 
solutions of the Schrodinger equation in the cylinder having boundary values 
of the form 

f{t,x)=a{t)h{x) 
where h G H"'-'^/'^{dMo). Given u G L2(M x Afo), we define 
lo 00 
Piv{t,x) =J2^it^^)'t>lix), P2v{t,x)= HtJ)Mx), 

1=1 l=lo+l 
where v{t,l) = • ), '/'^)L2{A/o)• 
Lemma 6.5. Assume that A G [Ai,oo) \ {Ai, A2, • • •}, choose /q ^ 1 so that 
Aio < A < A;„+i, let m > 2, and let Eq : H"'-^/'^{dMo) //'"(Mq) be a 
bounded extension operator. Let f{t,x) = a{t)h{x) where a G if™(M) with 
^ G M and G H'^-^/^idMo), and define 

Ef{t,x) = a{t)Eoh{x), Pf = -{-dj - Ag,+qo - A)F^. 
Also define 

Eji- ; f) = PjEf , Fj{--J) = PjPf. 
If ^ > 1/2, denote by u{t,x;f) the solution of the Dirichlet problem 

{-dt - Ago + go - A)n = in T, u\dT = f, 

{dt =F i\/A - \i)u{t, /;/)—;> as t -)• ±00 for ah 1 < / < /q, 
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where u{t,l]f) = {u{t, ■ ; f), cI}i)l2(^Mo)- Then we have 

U = Ul + U2, Uj = Ej + Wj, 

where Wj = Wj{- ; f) are the solutions of 

i-d^ - Ag, +qo- X)wi = ; /) in T, wi\eT = 0, 

{dt =F i\/ A — Xi)wi{t, /) — ;> as t ^ ±00 for ah 1 < / < Iq, 
with Wj{t,l) = {wj{t,-),(pi)L2i^Mo)^ and 

i-df - Ag, + qo - X)W2 = F2{- -J) inT, W2\aT = 0. 

If a G with 6 > 1/2, then the equation for wi has a unique solution 

wi G H!^s{T) n H^s,oiT) with u;i e Ran(Pi). Similarly, if a G H^{R) 
for some ;U G M, then the equation for W2 has a unique solution 1(72 G 
HJ^{T) n H'^ QiT) with 11)2 G Ran(P2)- We have the norm estimates 

ll^ilU™(T) + ll^illH™-2(r) - ^ll«llH;j»(R)ll^ll//"i-i/2(aAfo)' ^ ^) 

and 

II^^i||h™^(t) + I|w^i||h™^(t) < C|klU,7(M)ll^ll//™-i/2(aAfo)' > ^/^' 
Proof. We note the estimate 
Consequently 

The same estimates are true for Ej and Fj, since the projections Pj commute 
with dt, Agy and with multiplication by {t)^. The result now follows from 
Lemma 16.41 and the standard reduction from the Dirichlet problem to an 
inhomogeneous problem with zero boundary values. □ 

Proof of Proposition \6.2l Assume that A G [Ai, oo)\{Ai, A2, ...}, and choose 
^0 > 1 so that Xig < X < Aif,+i. Fix G M so that A — A;^ is not a Dirichlet 
eigenvalue of —Ag^ + qq in Mq. Also fix an integer m > dim(T)/2 + 1. 
We will show that for any h G //'"^^/^(SMq), one has the pointwise limit 

e'^'^Z.M - k')h = Jim ^ [^l,,M)fR' ^R' 
where fn is the function on dT given by 

fR{t,x) = e"''^R{t)h{x) 

and ^R{t) G C^(R) are suitable cutoffs. Below, we will use the notations 
in Lemma 16.51 For later purposes we choose the extension operator Eq : 
H-^~^/^{dMo) H'^iMo) in Lemma [631 to he Eo : h ^ v{x;h), where 
v{x;h) is the unique solution of the problem 

{-Agg+qo- {X- k'^))v{x;h) =0 in Mq, v{- ; h)\dMo = h. 

The proof below will make use of the splitting 

U = Ui+ U2, 
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and also the splitting 

V = Vl+V2 

where Vj{- ; h) G H'^{Mq) are the projections Vj{- ; h) = Qjv{- ; h). Here, for 
V G L2(Mo), 

lo oo 

Qiv = j2^mi{^), Q2V= ^ nO'AK^), 

«=1 l=lo+l 
with = (^^, </';)L2(Afo)- ^act, we will prove that 

lim -i- rd,uii-;fR>)dR' = e"'%vii-;h)\9T, (6.1) 

if,-S>oo K — I Ji 

hm -1- r d^U2i-;fR')dR' = e"'%V2i-;h)\aT. (6.2) 

iJ-5-oo K — I Ji 



Note that 



Thus the proposition will follow immediately from (|6.ip and (|6.2p . 

Let us next describe the cutoff functions. If -R — )• oo, the boundary value 
e*'^*^ij(t)/i(x) converges to e^^*h{x), a function in H]^ ^^"^{dT) for fi < —1/2. 
Fix some /i < -1/2, and let ipR{t) = 1 for \t\ < R, ipR{t) = for \t\ > R. 
We approximate the functions by ^R{t) e C^((-i? - 1, ii + 1)) that 

are functions for which 

lim W'^R - ijR\\Li(R) = 0, lim - 1||h™(m) = 0. (6.3) 

Such functions can be chosen to be 

1 for \t\ < R, 

$(i?"(|t| -i?)) iov \t\>R, 



^R{t) 



where $ G C^((— 1, 1)) is equal to one near 0, and a is a positive constant 
chosen so that ma + + 1/2 < 0. The norm bounds as i? — )• oo are valid 
because "^r — iI^r is supported in i? < |t| < i? + ^ r — 1 is supported 
in \t\ > R, W^rWw^,^ < CR""^, and {J^ t^^" dtY''^ < CR''+^/^. 
Let us denote 

fR{t, x) = e'^'-^R{t)h{x), fit, x) = e'^'h{x) 

where h G //'"^^/^(SMq). We will now prove (|6.2|) . Note that by construc- 
tion we have 

and thus Ej{- , Jr) — Ej{ - , f) = Ej{ - , Jr — f). The function Fj, wj, and Uj 
satisfy a similar property. Using Lemma 16.51 "we see that 

\\'W2i-;fR) - ii)2(- ;/)lk™(T) = \\'W2{-;fR - /)||H™(r) 

< C||e*^*(^'iJ - l)|U-(R)||/l||iym-l/2(5j^^g). 
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Similarly, 

\\E2{- ; Ir) - E2{- ; /)||//™(T) = \\E2{- ; /r - /)||h™(t) 

Since U2 = E2 + W2, the estimate ()6.3p implies that 
\\u2{- ; /r) - U2{- ; /)||//m(T) 

and consequently 

lim d^U2{- ; Jr) = duU2{- ; /) 

in ^^^(5T). By Sobolev embedding this limit also holds pointwise, and 
we also have 

1 

lim — — / d,U2{- ; /i?0 dB' = d,U2{- ; /)• (6.4) 

We can connect the last expression to the V2 component of the time- 
harmonic solution V = v{- ;h). As in Proposition l5.51 we see that the function 
e'-^^V2{x] h) G Hf^{T) solves the equation (-9| - A^^ + qo - \)u = in T 
with boundary value e^^^V2\dT- Similarly, the function U2 = U2{- ; /) solves 
the same equation with boundary value e*'^*((52-£'0^)|9T5 where Eq was the 
bounded extension operator. But since we chose EQh{x) = v{x;h), the two 
solutions have the same boundary values and by uniqueness one has 

U2{t,x-J) = e'^%2{x]h). 



Together with (|6.4p . this proves (|6.2p . 

It remains to show the identity ()6.ip concerning ui and vi. Recall that 



h 

ui{t,x;fii) = ^u{t,l;fR}(pi{x). 
1=1 



We start by giving formulae for Fourier coefHcients u{t, I; /r) for 1 < / < Iq. 
If 7/ £ H^^'^{dT) is compactly supported in the t variable, we have 

= {{-dt - Ago +qo- X)u{t, ■ ; 'n),(t>l)L^Mo) 



= {-dt + Xi-X)u{t,l;ri)+ / r^{t,y)d,My)dS{y). 

J a Mo 

Since A — A/ > for 1 < / < /q) this and the radiation condition imply that 

/oo 
Gi{t-t')fi{t\l)dt\ 
-00 

where 

GKt) = ;^(A-AO-^/Vl*l^, 

fi{t,l) = I r]{t,y)d,Uy)dSiy). 
J a Mo 
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1 /2 

Consider the function Ii{R) = Ii{t,x;R) in (dT) for any 6 > 1/2, 
given by 

lo 

Ii{R) = duUi{t,x;fR) = ^u{t,l; fii)du(t)i{x). 

1=1 

Using the expression for Fourier coefficients above, we have 



(=1 

^ (A - A,)-V2 



du4'i{x) 



E 

1=1 



2i 

(A - Az)-V2 



2i 



dT 



dT 



fR{t',y)d,Uy)dt'dS{y) 



^\'-''\'^'e"^''^R{t')h{y)d,Uy)dt'dS{y) 



du(f)i{x). 



Replacing here ^r by ipR and using the first estimate in (|6.3p results in an 
o(l) error in L^{dT) as R ^ oo. We thus obtain 



hiR) = Yl 



1=1 



2i 



-R 



h{i)d^(t>i{x) + o{i) 



where 



h{l)= / h{y)d,Uy)dS{y). 



Mo 



For a given t, we assume R so large that t G {—R,R). The t' integral can 
be computed explicitly, and we obtain 

^ i\t-t'\vy^ ikt' _ 2i(A - Xi^^^ ikt 



R 



J{k+^ \-Xi)R-iW p-iik-^\-Xi)R+iW\-\i 



+ 



i{k + VA^) - VX^i) 

The last two terms oscillate with respect to R, but we can remove these 
oscillating terms by averaging: since by assumption k it — A; 7^ 0, we 
have 



lim 



1 



R 



R^oo R — 1 

This shows that for any fixed (t, x) € (9T, we have 



1 1 

R^ooR-lJ^ ^ ^^X-k 



k^ 



We can relate the last expression to the time-harmonic solutions v{- ;h) of 
(— + go — (A — A;^))t> = in Mq with f laA/g = /i. We have 

= ((-^90 +90- {X-k^))v{-]h),(l)i)L2^Mo) 

= {\i - (A - A;2))t)(/; h)+ f h{y)d,(t>i{y) dS{y) 

J a Mo 
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which imphes 

h{l) = {\-\i-k^)v{l;h). 

This shows that 

•^^ 1=1 

= e''''d,vi{--h). 

This shows ()6.ip . which concludes the proof. □ 
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